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ABSTRACT 


The  chromatic  number  problem  is  to  determine  the  minimum 
number  of  colors  to  assign  to  the  vertices  of  a  graph  such 
that  no  connected  vertices  are  assigned  the  same  color. 
This  paper  presents  a  branch  and  bound  algorithm  for  the 
solution  to  the  chromatic  number  problem  and  investigates 
five  different  branching  functions.  Additionally#  a  method 
of  coloring  very  sparse  graphs  is  presented  which  divides  a 
graph  into  biconnected  components  and  reduces  the  time  re¬ 
quired  to  color  the  graph. 


3 


TABLE  OF  CONTENTS 


X.  INTRODUCTION  .  7 

II.  BACKGROUND  .  II 

A.  GRAPH  THEORY  .  11 

B.  BRANCH  AND  BOUND  .  14 

C.  PROBLEM  COMPLEXITY  .  19 

III.  THE  ALGORITHMS  . 27 

A.  GRAPH  GENERATION  .  27 

B.  GRAPH  DIVIDING  .  32 

C.  GRAPH  COLORING  .  38 

D.  THE  VERTEX  ORDERING  FUNCTIONS  .  45 

1.  Random . 47 

2.  Degree  . .  47 

3.  Dvec  3  . 4 

4.  Adjacency  . .  49 

5.  Least  Colors  Available  .  50 

E.  CONTROLLING  ALGORITHM .  52 

F.  GENERAL  COMMENTS  ON  ALGORITHM  DEVELOPMENT 

AND  THE  'C*  PROGRAMMING  LANGUAGE  .  54 

IV.  RESULTS  .  56 

A.  PARAMETERS  .  56 

1.  Time  . . 57 

2.  Space  . . 63 

B.  TIME  ANALYSIS . .  64 


1 


4 


C.  SPACE  ANALYSIS  .  68 

0.  ADVANTAGE  OF  OIVIDING  THE  GRAPHS  .  72 

V.  CONCLUSIONS  .  77 

APPENDIX  A— PROGRAM  SOURCE  CODE  .  80 

APPENDIX  B— TABULATION  OF  RESULTS  .  102 

LIST  OF  REFERENCES  .  122 

INITIAL  DISTRIBUTION  LIST .  12a 


I 


r  ■ 

* 


5 


LIST  OF  FIGURES 


1.  The  Tree  Structure  for  the  Colorability  Problem.  .  .  17 

2.  Comparison  of  Several  Polynomial  and 

Exponential  Complexity  Functions  ...  .  ..21 

3.  Transf  ormat  i  on  of  Problem  Q  to  Problem  R . 23 

4.  Examples  of  the  Isomorphs  of  2  Distinct  Graphs  ...  29 

5.  A  Graph  with  Articulation  Points . 32 

6.  A  Graph  Divided  at  its  Articulation  Points  .  33 

7.  A  Depth  First  search  of  a  Graph . 35 

8.  A  Schematic  of  a  Graph  with  3  Biconnected 

Components  and  2  Articulation  Points  .  36 

9.  Two  Paths  P  and  Q  which  are  not  Isomorphic  .....  44 

10.  Example  of  Shortcoming  of  Degree  Method  .  48 

11.  The  First  Four  Levels  of  the  Unpruned  Search  Tree.  .  60 

12.  Average  Expansions  vs.  Graph  Size  .  ....65 

13.  Maximum  Expansions  vs.  Graph  Size . 66 

14.  Average  Modes  vs.  Graph  Size.  .  69 

15.  Maximum  Nodes  vs.  Graoh  Size  . . 70 

16.  Average  Branching  Factor  vs.  Depth  in  the  Tree  ...  71 

17.  Ratio  of  Largest  Biconnected  Component 

to  the  Original  Graph  Size  . . . . 73 

18.  Ratio  of  Expansions  Required  for  Divided  Graph 

to  Expansions  Required  for  Undivided  Graph  .  75 


6 


I.  INTRODUCTION 


Informally*  a  graph  is  a  collection  of  points#  called 
nodes  or  vertices*  some  of  which  may  be  connected  by  lines* 
called  edges.  A  solution  to  the  graph  colorability  problem 
specifies  a  color  for  each  node  such  that  any  two  nodes 
which  are  connected  by  an  edge  do  not  have  the  same  color. 
The  optimal  colorability  problem  is  to  assign  colors*  as  in 
the  graph  colorability  problem#  with  the  restriction  that 
the  fewest  number  of  colors  possible  is  used. 

For  centuries*  mathematicians  have  studied  graphs  and 
their  characteristics  with  the  awareness  that  many  engineer¬ 
ing  and  mathematical  problems  could  be  modeled  by  a  graph. 
Then  questions  about  the  engineering  problem  could  be 
answered  by  studying  the  characteri st ics  of  the  graph  and 
using  knowledge  derived  from  graph  theory. 

One  of  the  earliest  and  best  documented  graph  colorabil¬ 
ity  problems  dates  back  to  1852#  when  a  student  of  De  Morgan 
named  Francis  Guthrie#  proposed  that  4  colors  was  all  that 
was  needed  to  color  a  map  such  that  no  two  bordering  coun¬ 
tries  are  the  same  colorIRef.  11.  Note  that  this  problem  is 
transposed  into  a  graph  colorability  problem  by  using  a  node 
for  each  country*  and  connecting  2  nodes  with  an  edge  if  the 
two  correspondi ng  countries  share  a  common  border. 


1 

1 


7 


In  fact/  beginning  in  the  19th  century#  many  of  the  en¬ 
gineering  problems  in  areas  such  as  scheduling#  communica¬ 
tions#  transportat ion#  electronics#  and  chemistry  have  been 


studied  through  the  use  of  graph  theory.  More  recently#  the 
social#  biological#  and  environmental  sciences  have  relied 
on  graph  theory  to  investigate  specific  problems.  In  the 
decade  of  the  seventies#  computer  science  and  graph  theory 
have  been  closely  intertwined  as  computer  scientists  real¬ 
ized  that  many  of  their  difficult  proolems  can  best  be  un¬ 
derstood  by  graph  modeling.  One  of  the  most  notable  is  the 
scheduling  problem  where  two  or  more  processes  must  compete 
for  common  resources.  Ways  to  efficiently  schedule  these 
processes  are  better  understood  when  the  system  is  modeled 
by  a  graph. 

The  intent  of  this  research  has  been  to  take  a  single 
Droblem  from  graph  theory#  namely  optimal  graph  colorabili- 
ty»  and  to  investigate  the  use  of  the  digital  computer  to 
find  optimal  colorings  for  various  graphs.  Though  a  seem¬ 
ingly  trivial  problem  on  a  computer#  it  is  complicated  by 
the  fact  that  in  the  worst  case#  all  known  methods  of  find¬ 
ing  solutions  to  problems  of  this  type  reouire  exponential 
time.  That  is#  as  the  size  of  the  graph  gets  larger#  the 
time  required  to  find  the  optimal  coloring  grows  exponen¬ 
tially. 
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Typically#  the  researcher  investigating  an  exponential 
time  problem#  whose  size  is  very  large#  is  faced  with  one  of 
two  choices:  1)  he  must  accept  a  solution  which  is  not 
necessarily  optimal#  but  some  approximate  or  'close*  solu¬ 
tion#  or  2)  he  must  relax  some  of  the  restrictions  on  the 
problem.  Most  of  the  research  done  on  the  graph  colorabili- 
ty  problem  ias  been  in  the  area  of  finding  approximation  al¬ 
gorithms#  that  is#  algorithms  which  attempt  to  find  a  rela¬ 
tively  few  number  of  colors  without  undue  computational  ef¬ 
fort  being  expended  [Ref  a .  2»3#4],  The  primary  thrust  of 
this  research  has  been  to  investigate  methods  of  finding 
only  optimal  colorings#  and  to  study  the  effects  of  dif¬ 
ferent  procedures  on  the  time  and  space  (computer  memory) 
requi red  by  each  . 

Specifically#  the  branch  and  bound  technique  is  used  to 
find  the  ootimal  coloring  of  a  large  number  of  randomly  gen¬ 
erated  graphs#  and  the  branching  function  is  varied  to 
determine  the  relative  value  of  different  heuristic  Branch’ 
i ng  functions. 

In  summary#  this  research  encompasses  three  related#  yet 
diversified  disciplines.  From  mathematics  comes  graph 
theory#  from  operations  research  comes  the  branch  and  bound 
techniques;  and  from  computer  science  comes  problem  complex¬ 
ity.  Background  information  is  provided  for  these  three  to- 
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pics  in  the  first  Dart  of  this  report.  Following  that  is  a 
discussion  of  the  research  done  on  the  computer#  including 
the  algorithms  used#  the  branching  functions  investigated# 
and  the  graph  generation  techniques.  Finally#  there  is  a 
discussion  of  the  results  obtained#  the  difficulties  encoun¬ 
tered#  and  the  ideas  and  notions  derived  during  the  course 
of  the  research. 
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II .  8ACKGR0UND 


A.  GRAPH  THEORY 

While  not  attempting  in  any  way  to  fully  cover  the  field 
of  graph  theory*  an  introduction  to  the  notions  and  termi¬ 
nology  which  relate  to  the  coloraDility  problem  is  necessary 
for  an  understanding  of  the  discussion  which  follows.  The 
definitions  used  here  are  standard  in  most  of  the  literature 
on  graph  theory*  so  the  reader  with  some  background  in  the 
field  can  skip  this  section  without  any  loss  of  continuity. 

Formally*  a  graph  G  =  (V*£)  consists  of  a  finite* 
nonempty  set  of  nodes  or  vertices  V*  and  a  finite*  possibly 
empty  set  of  vertex  pairs  called  edges  E.  For  the  purposes 
of  this  research  the  set  of  edges  is  restricted  to  not  con¬ 
tain  any  elements  of  the  form  (x*x).  That  is*  no  edge  can 
exist  from  a  vertex  to  itself.  In  the  case  where  the  edges 
have  direction*  they  are  defined  by  an  ordered  pair  of  ver¬ 
tices  and  the  graph  is  called  a  DIGRAPH.  This  research  has 
been  strictly  limited  to  non-directed  graphs. 

The  general  graph  colorability  problem  is  formalized  as 
follows:  Given  graph  G  and  integer  k*  does  there  exist  a 
function  R:  V  ->  I  where  I  =  <1*2*. ..*k>  such  that  for  all 
(x,y)  S  E*  R(x)  t  R(y)?  The  optimal  colorability  problem  is 
to  find  the  smallest  integer  k  such  that  a  function  R  ex¬ 
ists.  That  smallest  integer  found  is  called  the  chromatic 


number  of  t he  graph  G.  If  R(x)  -  R(y)  then  vertices  *  and  y 
are  said  to  belong  to  the  same  color  class.  Note  that  the 
function  R  partitions  the  vertices  of  the  graph  into  k  color 
classes. 

Some  additional  terms  relating  to  graphs  are: 

Adjacent:  Vertex  x  is  adjacent  to  vertex  y  if  an  edge  ex¬ 

ists  between  them#  that  is  (x#y)  6  E. 

Adjacency  matrix:  A  two  dimensional#  n  x  n  matrix#  used  to 
describe  any  graph.  The  indices  into  the  matrix  are 
the  vertices  of  the  graph#  and  the  elements  are  a  1 
if  an  edge  exists  between  the  two  vertices#  and  a  0 
otherwise.  For  the  graphs  used  in  this  research  the 
diagonal  elements  are  all  zeroes  and  the  matrix  is 
reflective  about  the  diagonal.  The  adjacency  matrix 
is  a  very  common#  and  easy  to  implement  method  of 
representing  a  graph  in  a  computer. 

Degree:  The  degree  of  a  vertex  is  the  number  of  edges  in¬ 

cident  upon  it#  or  the  number  of  vertices  which  are 
adjacent  to  it.  The  degree  of  vertex  x  is  easily 
calculated  from  the  adjacency  matrix  by  simply  sum¬ 
ming  the  l's  in  row  x  or  column  x.  If  the  elements 
of  the  adjacency  matrix  are  denoted  a(i#j)  and  the 
degree  of  vertex  x  is  denoted  d(x)  then 

d(x)  =  X  a  ( x  #  i  ) 
i 
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Path:  A  oath  is  said  to  exist  from  vertex  x  to  another  ver¬ 

tex  y  if  one  can  start  at  vertex  x  and  folio*#  the 
edges  of  the  graph  and  reach  vertex  y.  Formally,  a 
path  exists  if  there  exists  an  ordering  of  some  of 
the  vertices  of  G  ( V} » V2 , . .  •  , vn  )  such  that  v^  =  x  * 
vR  =  y  and  for  1  <=  i  <=  n-l»  (v^/v^j)  6  E. 

Connected  Graoh:  A  graph  is  connected  if  at  least  one  path 
exists  from  each  vertex  to  every  other  vertex  in  V. 

Doubly  connected  graoh:  A  graoh  is  doubly  connected  if  for 
every  vertex  x  and  every  vertex  y  at  least  2  oaths 
exist  from  x  to  y  and  no  other  vertices  are  common 
to  both  oaths.  This  is  also  called  a  biconnected 
graph . 

Articulation  point:  A  vertex  v  is  an  articulation  point  of 
graph  G  if  there  exists  some  vertex  x  and  some  ver¬ 
tex  y  such  that  every  oath  from  x  to  y  includes  ver¬ 
tex  v.  Note  that  a  doubly  connected  graph  has  no  ar¬ 
ticulation  points. 

Complete  graph:  A  graph  G  is  complete  if  all  possible  edges 
are  present.  That  is  G  is  complete  if  for  all  ver¬ 
tices  x  and  y,  (x,y)  6  E. 

Cligue:  A  clique  is  a  subset  of  a  graph  such  that  every 

vertex  in  the  clique  is  adjacent  to  every  other  ver¬ 
tex  in  the  c 1 i que . 
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B.  BRANCH  AND  BOUND 

Many  problems  which  deal  with  searching  for  a  solution 
or  set  of  solutions  satisfying  some  constraints  can  be 
solved  using  the  branch  and  bound  technique.  It  is  espe¬ 
cially  useful  in  solving  minimization  (or  maximization) 
problems  and  has  been  successfully  used  in  constrained  op¬ 
timization  Drool ems  since  the  late  1950's.  Although  there 
are  several  ways  of  describing  branch  and  bounds  the  basic 
idea  is  to  split  the  solution  space  (branch)  and  place  a 
limit*  or  lower  bound*  on  the  optimal  cost  of  the  problem 
limited  to  each  subset  of  the  solution  space.  Those  subsets 
whose  optimal  costs  do  not  exceed  the  cost  of  some  known 
solution  are  then  repeatedly  divided  ana  bounded  until  a 
solution  is  found  whose  actual  cost  is  no  greater  than  the 
lower  bound  on  any  of  the  subsets.  This  solution  is  thus 
the  optimal  one.  The  power  of  branch  and  bound  comes  from 
its  ability  to  leave  unexplored  those  suosets  which  are 
known  not  to  contain  the  optimal  solution. 

In  order  to  use  the  branch  and  bound  technique  the  prob¬ 
lem  must  have  an  expressable  set  of  solutions  and  a  cost 
function  COST()  which  maps  the  solutions  into  nonnegative 
inteqers.  Usually*  the  problem  calls  for  finding  that  solu¬ 
tion  s  for  which  COST(s)  is  a  minimum.  Sometimes  all  solu¬ 
tions  may  be  desired  whose  cost  is  no  greater  than  COST(s). 


Many  problems  require  that  all  the  solutions  satisfy 


some  set  of  constraints  which  may  be  divided  into  two 
categories*  explicit  and  implicit.  If  a  solution  is 
represented  by  an  n-tuple  ( x^  *  •  •  •  •  •  *n  )  *  then  the  explicit 

constraints  are  ones  which  determine  what  values  the  x..  's 
may  take  on.  The  implicit  constraints  restrict  the  ways  in 
which  the  x. 's  in  a  solution  relate  to  each  other.  All  solu¬ 
tions  which  satisfy  the  explicit  constraints  are  said  to  oe- 
long  to  the  possible  solution  space.  The  implicit  con¬ 
straints  then  determine  which  solutions  in  the  solution 

space  satisfy  the  constraints  of  the  problem.  These  solu¬ 
tions  make  up  the  'feasible'  solution  space. 

For  the  graph  colorability  problem*  we  will  use  the  con¬ 
vention  that  a  solution  is  expressed  as  an  n-tuple 

( x,  * x  *  .  .  .  * x  )  of  integers  where  l<  =  x.<sn*  and  x.  is  the 

1  2  n  l  l 

color  assigned  to  vertex  i.  The  explicit  constraints  res¬ 
trict  the  x. 's  to  integer  values  from  1  to  n.  Therefore* 
the  possible  solution  space  has  nn  elementsCn  ways  of  pick¬ 
ing  an  integer  from  1  to  n  for  each  of  n  positions).  The 
implicit  constraint  in  this  problem  is: 

(i* j)  6  E  =>  x. *  x. 

I  J 

The  size  of  the  feasible  solution  space  is  therefore 
also  a  function  of  the  number  of  edges  in  the  graph  and 


their  arrangement 


The  first  steo  of  branch  and  bound  is  to  systematically 


divide  the 

sol ut i on 

space  into 

subsets . 

This 

may 

be 

represented 

by  a  tree 

organ i zat i on , 

where  the 

root 

of 

the 

tree  reoresents  all  feasible  solutions.  The  children  of  a 
node  N  then  reoresent  the  subsets  into  which  N  can  be  divid¬ 
ed.  The  leaf/  or  terminal  nodes  are  the  particular  solu¬ 
tions  which  make  up  the  set  of  feasible  solutions. 

For  the  graph  colorability  problem,  the  following  tree 
organization  was  used.  The  root  node  represents  all  feasi¬ 
ble  n-tuples.  Then,  an  arbitrary  vertex  is  chosen  and  the 
n-tuoles  are  divided  into  n  subsets  such  that  subset(k)  con¬ 
tains  all  the  n-tuples  for  which  x^  s  k.  This  is  then  done 
repeatedly,  picking  a  new  vertex  at  each  level  of  the  tree. 
For  notation  purposes,  a  partial  solution  ( x^  ,  x^  ,  • . . ,  x^ )  , 
k  <  n,  is  an  assignment  of  colors  to  the  vertices  1,2,  ...,k. 
This  k-tuple  is  then  used  to  reoresent  all  the  elements  of 
a  subset.  Using  this  notation,  a  portion  of  the  tree  for 
the  case  n  =  4  is  shown  in  fig.  1.  Note  that  the  lower 
bound  on  the  cost  of  a  node  is  the  number  of  different 
colors  in  the  partial  solution. 

A  branching  strategy  is  a  rule  for  determining  how  a 
solution  is  to  be  divided  into  subsets.  For  this  problem 
the  branching  function  is  the  method  of  choosing  the  next 
vertex  to  color  at  each  level  of  the  tree. 


(  ) 


(1,2,2, 1)  a ,2,2,2; 


The  tree  structure  for  the  Colorability  Problem. 

Figure  i 


Once  the  organisation  is  decided  upon,  the  task  of  the 
problem  solving  method  is  to  explore,  or  search  this  'state 
space  tree'  until  the  optimal  answer  state  is  found.  Since 
the  time  required  by  the  search  algorithm  is  a  function  of 
the  number  of  nodes  explored  in  the  tree,  the  search  stra¬ 
tegy  is  to  explore  the  fewest  nodes  possible.  One  way  to 
facilitate  this  is  to  'prune'  the  tree  as  early  in  the 
search  as  Possible.  Pruning  takes  place  when  a  node  is 


reached  which  violates  the  implicit  constraints  of  the  orob 


lem.  Since  no  legal  answer  state  could  oossibly  be  reached 
through  that  node*  none  of  its  children  need  be  explored. 
The  size  of  the  search  soaee  is  thus  reduced  by  the  number 
of  nodes  below  the  pruned  node.  Obviously*  the  higher  in 
the  tree  that  pruning  takes  place*  the  greater  the  reduction 
in  the  search  soace. 

The  second  way  to  reduce  the  search  space  is  the  essence 
of  the  branch  and  bound  method.  Each  node  has  a  cost  asso¬ 
ciated  with  it*  and  since  the  cost  function  is  non¬ 
decreasing*  no  searching  need  take  place  below  any  node 
whose  cost  is  greater  than  the  cost  of  some  known  solution. 
A  'search  strategy'  is  a  rule  for  choosing  which  unexolored 
node  to  explore  next  in  the  tree  search.  Several  different 
basic  search  strategies  exist  for  exploring  a  tree. 

A  deoth  first  search  proceeds  down  the  tree  from  the 
root*  exploring  the  leftmost  children  of  each  node  first. 
At  any  point  where  pruning  takes  place*  the  search  is  backed 
uo  to  the  first  node  where  it  can  again  proceed  down  the 
tree  in  a  different  oath.  Too-to-bottom*  1 ef t-to-r i ght  is  a 
good  description  of  this  method. 

A  breadth  first  search  explores  all  the  children  of  a 
node  before  oroceedina  down  the  tree.  This  might  be  called 
a  1 ef t-to-riqht *  top-to-bot tom  search. 


A  'best-first'  or  'least  cost'  search  strateqy  e*o)ores 
only  the  most  promising  path  throuqh  the  tree  at  any  qiven 
time.  A  node  is  explored  by  assigning  costs  to  all  of  its 
children  and  placing  them  on  an  unexplored  list.  The  next 
node  to  be  explored  is  the  one  on  the  unexplored  list  whose 
cost  is  currently  the  lowest.  A  priority  queue  can  be  used 
to  maintain  an  orderinq  of  the  unexplored  nodes  and  at  the 
point  where  the  highest  priority  node  is  an  answer  state 
then  the  search  can  be  terminated  because  that  node 
represents  the  best  or  optimal  solution  to  the  problem. 
Nodes  which  violate  the  problem  constraints  may  be  assigned 
an  arbitrarily  high  cost/  or  merely  not  olaced  in  the  queue# 
which  effectively  prunes  the  tree.  This  is  the  search  stra¬ 
tegy  employed  in  this  research  effort. 

This  has  been  a  general  description  of  the  branch  and 
bound  technique.  More  specific  details  of  the  tree  arrange¬ 
ment#  the  search  strategy#  and  the  cost  function  used  in 
this  research  will  be  described  in  Section  III. 

C.  PROBLEM  COMPLEXITY 

The  complexity  of  a  problem  is  said  to  be  polynomial- 
time  if  an  algorithm  exists  for  which  the  number  of  elemen¬ 
tary  or  fundamental  operations  is  limited  by  a  polynomial  of 
the  length  of  an  encoding  of  the  problem.  A  problem  is 
oo I ynomi a  I -space  if  an  algorithm  exists  for  which  the  amount 


of  elementary  computer  storage  space  needed  for  computation 
is  never  greater  than  some  polynomial  function  of  the  encod¬ 
ing  of  the  problem.  An  algorithm  for  a  ool ynomi al -t i me 
problem  is  said  to  be  a  'good*  algorithm  and  all  problems 
for  which  a  good  algorithm  exist  belong  together  in  a  class 
called  P.  It  may  seem  somewhat  extravagant  to  group  all 
ool ynomi al -t i me  algorithms  into  one  class  since  polynomials 
can  be  quite  large.  However/  no  matter  what  the  coeffi¬ 
cients  are/  in  the  limit  as  n  gets  large/  every  exponential 
function  dominates  every  polynomial  function.  Further/ 
though  unexplained/  experience  has  shown  that  for  many  of 
the  problems  encountered  which  have  a  oo 1 ynom i a  1 -t i me  algo¬ 
rithm/  the  solution  is  bounded  by  a  polynomial  of  small  de¬ 
gree.  Some  examples  are  ordered  searching  which  is  0(ln  n), 
sorting  which  is  0(n*ln(n))/  and  matrix  multiplication/ 
which  is  0(n^ ' ) . 

Another  group  of  problems  are  those  whose  best  'known' 

algorithm  require  greater  than  polynomial  time.  Some  of 

these  are!  1)  Algorithms  requiring  subexponential  time/  such 

as  0(e'^'")z  2)  algorithms  requiring  exponential  time/  such  as 

0(en)z  and  3)  algorithms  requiring  super-exponent i a  1  time/ 
en 

such  as  0(e  )•  The  disadvantage  of  non-po I ynomi a  1 -t i me  al¬ 

gorithms  is  in  the  explosive  growth  of  the  maximum  solution 
time  as  illustrated  by  figure  2.  In  this  table  the  maximum 
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computing  time  is  shown  for  different  complexity  functions 
as  n  goes  from  10  to  60. 


F  i  me 

compl ex i t  y 
funct i on 

10 

20 

Size 

30 

n 

40 

50 

60 

n 

.00001 

second 

.00002 

second 

.00003 

second 

.00004 

second 

.00005 

second 

.00006 

second 

n2 

.0001 

second 

.0004 

second 

.0009 

second 

.0016 

second 

.0025 

second 

.  0036 
second 

n3 

.001 

second 

.008 

second 

.027 

second 

.064 

second 

.125 

second 

.216 

second 

n5 

.  1 

second 

3.2 

seconds 

24.3 

seconds 

1.7 

mi nutes 

5.2 

mi  nutes 

13.0 

minutes 

2n 

.001 

second 

1.0 

second 

17.9 
mi nutes 

12.7 

days 

35.7 

years 

366 

centuries 

3n 

.059 

second 

58 

mi nutes 

6.5 

years 

3855 
cent . 

2xl08 

cent . 

1.3xl0l3 
cent  . 

Comparison  of  several  polynomial  and  exponential  time 
complexity  f unct i ons  [Ref .51 
Figure  2 


Note  that  the  only  concern  is  with  known  algorithms  and 
worst-case  situations.  No  conclusions  are  made  about  ex¬ 
pected  or  average  performance  of  an  algorithm. 

Problems  for  which  there  exists  an  algorithm  which  can 
guess  a  solution  and  verify  its  correctness  in  polynomial- 
time  make  up  the  class  of  non-det erm i n i s t i c  pol ynom i a  1 -t i me 
problems*  NP.  For  example*  graph  coloring  is  in  NP  because 
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an  algorithm  exists  which  can  guess  a  coloring  and  check 
whether  or  not  it  is  legal  in  polynomial-time.  Such  an  al¬ 
gorithm  needs  0(n)  time  to  arbitrarily  assign  a  color  to 
each  of  the  n  vertices.  Then  it  needs  to  consider  each  edge 
of  the  graph  and  check  whether  the  2  vertices  have  the  same 
color.  This  can  be  done  in  constant  time  for  each  edge  and 
since  there  may  be  0(n2)  edges*  the  entire  verification 
takes  0(n2) . 

Note  in  all  of  this  discussion  the  word  'known*.  There 
is  not  one  single  problem  in  the  class  NP  for  which  it  has 
been  proved  that  the  oroblem  is  not  in  Ptftef.  6],  There  are 
no  po l ynomi a  1 -bounded  algorithms  known  for  many  problems  in 
NP,  but  no  lower  bound  which  is  greater  than  polynomial  time 
has  been  proven  for  any  of  these  problems.  Thus*  it  is 
still  an  open  guestion  as  to  whether  the  class  NP  properly 
contains  P  or  whether  the  two  classes  are  equal. 

There  is  obviously  much  interest  in  trying  to  prove 
whether  a  qiven  problem  Q  could  or  could  not  be  solved  in 
ool ynomi al -t i me.  For  if  one  can  prove  that  no  polynomial¬ 
time  algorithm  for  Q  could  possibly  exist  then  there  is  no 
point  in  expending  the  effort  to  find  such  an  algorithm. 
Short  of  proving  that  a  oroblem  has  no  oo l y nom i a  1 -t i me  algo¬ 
rithm*  however*  there  is  some  comfort  in  knowing  that  ones 
difficult  oroblem  is  somehow  related  or  at  least  as  diffi- 
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cult  as  many  other  problems  for  which  no  one  has  found  a 
Pol ynomi al -t i me  algorithm.  For  this  discussion  some  more 
definitions  are  needed.  A  'problem'  is  some  general  ques¬ 
tion  to  be  answered  which  usually  has  one  or  more  unspeci¬ 
fied  parameters  and  some  specified  properties  which  the 
answer  is  required  to  satisfy.  An  'instance'  of  a  problem 
is  obtained  when  the  parameters  are  specified.  For  example, 
an  instance  of  the  colorability  problem  is  obtained  by 
specifying  the  vertices  and  edges  of  the  graph  and  the 
number  of  colors  for  which  a  coloring  is  desired. 

A  problem  Q  is  said  to  be  'reducible'  or  'transform¬ 
able*  to  R  if  an  instance  of  0  can  be  transformed  in  poly¬ 
nomial  time  to  an  instance  of  R*  In  other  words.  if  the 
answer  to  problem  Q  is  'yes'  if  and  only  if  the  answer  to  R 
is  'yes'  and  if  the  input  parameters  to  problem  R  can  be 
determined  in  po 1 ynom i a  1 -t i me  given  the  input  parameters  to 
Q»  then  Q  is  reducible  to  R.  The  following  schematic  illus¬ 
trates  this  definition: 


Input 
to  0 


Po 1 ynom i a  1 

Input 

A l gor i t  hm 

t  i  me 

to  R 

t  ransf orm 

for  R 

to  R 

Output 
for  3 
and  Q 


> 


Trans f ormat i on  of  Problem  Q  to  Problem  R 

Figure  3 
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Clearly*  if  Q  is  reducible  to  R  and  R  is  in  the  class  P* 
then  Q  is  also  in  P.  Note  that  nothing  is  said  about  the 
existence  of  an  algorithm  for  R*  only  that  'if'  an  algorithm 
for  R  exists  which  runs  in  doI ynomi al -t i me*  then  one  also 
exists  for  Q. 

If  every  problem  in  the  class  NP  is  reducible  to  some 
problem  Q  then  Q  is  said  to  be  NP-hard.  In  other  words*  Q 
is  at  least  as  hard  as  every  problem  in  NP.  Further*  if  Q 
is  also  in  the  class  NP  then  Q  is  NP-complete.  The  signifi¬ 
cance  of  a  oroblem  being  NP-complete  is  that  if  a 
po 1 ynom i a  1 -t i me  algorithm  is  ever  found  for  any  NP-complete 
oroblem  then  it  will  be  known  that  a  po 1 y nom i a  1 -t i me  algo¬ 
rithm  exists  for  all  NP  problems  and  that  indeea  NP  =  P. 
Further*  if  any  oroblem  in  NP  is  ever  shown  to  require 
suoer-ool ynomi al  time  complexity  then  every  NP-complete 
Droolem  must  also  reoui re  greater  than  polynomial  time. 

To  show  that  a  particular  problem  Q  is  NP-complete* 
then*  one  must  show  that  it  is  indeed  in  NP  ana  then  either 
show  that  every  NP  problem  could  be  reduced  to  Q  or  that 
some  already  proven  NP-comolete  problem  is  reducible  to  Q. 
In  1971*  Steven  Cook  (Ref.  7]  laid  the  ground  work  for  the 
current  theory  of  NP-comp 1 et eness  in  a  paper  entitled  "The 
Complexity  of  theorem  proving  procedures."  One  of  his  most 
significant  results  was  the  proof  that  the  'Satisfiability' 


oroblem  from  boolean  logic  is  NP-compl ete*  and  it  holds  the 
distinction  of  being  the  first  NP-comolete  oroblem.  Refer¬ 
ence  5  contains  an  excellent  description  of  the  Satisfiabil¬ 
ity  problem  as  well  as  a  proof  of  Cook's  theorem.  Cook  also 
proved  that  a  variation  of  the  Satisfiability  problem* 
called  3-Satisfiability*  was  NP-complete  by  transforming  Sa¬ 
tisfiability  into  3-Sat i s f i ab i 1 i t y  .  Then,  in  1972* 
KarplRef.  81  proved  that  many  decision  problems  were  also 
NP-complete.  One  of  these  proofs  was  that  the  3- 
Satisfiability  oroblem  was  reducible  to  the  colorability 
problem.  Thus  the  colorability  problem  is  also  NP-comolete. 

Inasmuch  as  most  theorists  are  in  agreement  that  the 
NP-complete  oroblems  are  probably  suoer-po 1 y nom i a  1  *  much 
work  has  been  done  in  recent  years  to  find  relaxations  to 
the  problems  which  allow  them  to  be  solved  in  polynomial 
time.  As  the  graph  colorability  proolem  has  many  apolica- 
tions  in  real-life  problems*  many  methods  of  finding  approx¬ 
imate  solutions  have  been  tried.  In  many  cases  where  a 
graph  can  be  used  to  model  a  real-life  problem*  the  absolute 
minimum  number  of  colors  may  not  be  needed*  but  rather*  some 
close  approx i mat i on  might  be  sat i sf actory  .  Unfortunately* 
even  this  has  proved  to  be  a  difficult  oroblem.  Many  algo¬ 
rithms  have  been  dev i sed  (Re f s  .  2* 3*u]  to  find  approximate 
colorings*  but  for  every  one*  it  is  possible  to  construct  a 
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graph  such  that  the  approximation  algorithm  will  find  a 
coloring  which  uses  many  more  colors  tfian  the  chromatic 
number.  Though  the  algorithm  may  indeed  run  in  polynomial 
time#  the  closeness  of  the  approximation  to  the  actual 
chromatic  number  of  the  graph  may  be  suspect.  In  fact#  in 
1976#  Garey  and  Johnson  IRef .  1^1  proved  the  following: 

Given  graph  G  with  chromatic  number  X(G)#  and 
do  1 ynom i a  1 -t i me  approximation  algorithm  A#  which  computes  a 
coloring  of  G  using  A(G)  colors#  there  can  be  no  A  such  that 
for  all  G : 


A  ( G ) 

"  "  <r  for  r  <  2, 

X  ( G 1 

In  other  words#  any  approximation  algorithm  will  always 
have  some  input  graph  for  which  it  can  find  no  coloring 
which  uses  less  than  twice  as  many  colors  as  actually  need¬ 
ed. 

The  primary  reference  for  the  material  in  this  section 
is  Garey  and  JohnsonlRef.  53.  It  is  highly  recommended  for 
the  reader  who  is  interested  in  other  NP-complete  problems 
and  further  discussion  on  this  topic. 
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III.  THE  ALGORITHMS 


A.  GRAPH  GENERATION 

One  of  the  earliest  problems  faced  in  this  research  Mas 
that  of  deciding  what  ch arac t er i s t i c s  Mere  desirable  in  the 
qraphs  to  be  colored  and  how  to  generate  them.  As  the  goal 
of  the  research  Mas  to  contrast  the  efficiencies  of  various 
heuristics  for  optimal  colorings#  the  graphs  used  needed  to 
be#  in  some  sense#  typical#  or  representative  of  graphs 
Mhich  occur  in  real-life  problems.  Though  the  question  of 
Mhat  is  a  tyDical  graph  Mas  left  largely  unanswered#  three 
possibilities  for  graph  generation  Mere  considered. 

One  method  considered  Mas  to  use  hand  generated  graphs 
Mith  certain  built-in  specific  features.  For  example# 
graphs  could  be  designed  Mith  particular  distributions  of 
degree  of  the  vertices.  Another  method  Mould  be  to  generate 
graphs  Mith  a  given  chromatic  number  but  Mith  varying 
numbers  of  vertices  or  edges.  The  major  draMbacks  to  this 
method  Mere  the  necessarily  small  sample  of  graphs  and  the 
failure  of  the  graphs  to  be  typical  in  any  sense.  This 
method  Mas  therefore  not  used  for  any  of  the  experiments 
done  but  Mas  used  extensively  for  program  testing  and  debug¬ 
ging. 

Another  possibility#  quickly  discarded#  Mas  to  survey 
some  of  the  disciplines  Mhere  graph  colorability  models  the 


problem  and  to  use  some  actual  graphs  which  occur  in  real- 
life.  Though  this  might  generate  some  interesting  results 
for  those  specific  applications#  it  would  also  lead  to 
results  which  were  heavily  flavored  by  the  problem  peing 
considered.  what  was  needed  was  a  selection  of  graphs  which 
were*  in  some  sense#  represent  at i ve  of  the  set  of  all 
graphs  . 

The  most  promising  method  of  graph  generation  was  to  use 
the  computer  and  to  create  a  large  population  of  random 
graphs  on  which  to  test  the  various  heuristics.  There  is 
some  benefit  to  using  random  graphs  since  one  generally 
feels  comfortable  that  these  random  graphs  are  fairly 
representative  of  the  set  of  all  graphs. 

The  primary  drawback  to  random  generation  of  graphs  is 
that  the  generator  will  tend  to  generate  more  of  the  types 
of  graph  which  have  many  isomorphisms  than  the  type  with  few 
isomorphisms.  A  graph  G  is  isomorphic  to  qraph  H  if  the 
vertex  numbers  of  G  can  be  permuted  in  such  a  way  as  to  make 
the  resulting  adjacency  matrix  of  G  identical  to  the  adja¬ 
cency  matrix  of  H.  Fig  4  shows  2  of  the  isomorphic  graphs 
for  each  of  the  only  2  distinct  graphs  possible  for  4  ver¬ 
tices  and  4  edges. 

Since  graph  coloring  algorithms  can  be  influenced  by  the 
different  types  of  graphs#  a  random  generator  may  bias  the 
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investigation  of  the  relative  performances  of  different 
coloring  programs.  The  generation  technique  could  be  amend¬ 
ed  by  somehow  reducing  the  probability  of  using  a  particular 
graph  by  dividing  the  probability  by  the  number  of  the 
graph's  isomorphisms.  But  determining  the  number  of  iso¬ 
morphic  configurations  of  a  given  graph  is  not  a  straight¬ 
forward  problem  and  beyond  the  scope  of  this  research. 


Examples  of  the  isomorphs  of  2  distinct  graphs. 

Figure  4 

Another  technique  might  be  to  record  each  graph  and  to 
test  each  new  one  generated  against  all  previous  to  deter¬ 
mine  whether  it  is  an  isomorphism  of  an  already  used  graph. 
But  the  problem  of  just  testing  two  graphs  against  each  oth¬ 
er  is  probably  NP-compl et e ( t hough  considered  an  open  ques¬ 
tion  at  this  time)  £Ref.  51. 

In  order  to  get  on  with  the  research  goals  at  hand»  it 
was  decided  to  accept  a  random  graph  generator  which  gen¬ 
erates  graph  type  G  more  than  type  H  if  G  has  more  isomor¬ 
phisms  than  H,  This  is  an  admitted  shortfall  in  this 
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research  effort  and  it  is  booed  that  the  algorithms  experi¬ 


mented  with  are  not  discriminated  against  by  the  shortcoming 
of  this  random  graph  generator. 

Inasmuch  as  unconnected  graphs  are  no  more  difficult  to 
color  than  their  largest  connected  component*  it  was  decided 
to  run  the  initial  tests  on  only  connected  graphs.  Two 
methods  of  generating  connected  graphs  were  considered:  1) 
To  randomly  add  edges  until  the  graph  was  connected  and  2) 
to  create  a  spanning  tree  of  n  vertices  and  n-1  edgesithe 
minimum  to  connect  the  n  vertices)  and  then  add  additional 
edges  at  random.  Since  method  one  would  require  testing  tor 
connectivity  after  each  edge  and  also  would  fail  to  generate 
very  many  sparse  graohstfound  often  in  real  oroblems)  method 
2  was  chosen  and  implimented  using  the  following  algorithm: 


Algorithm  GENGRAPH 

1.  Let  VC  be  a  set  of  connected  vertices*  initially  empty 
Let  VU  be  a  set  of  unconnected  vertices*  initially 

VU  -  (1*2*. . • * n ) 

2.  Choose  a  random  vertex  a 

3.  Move  a  from  VU  to  VC 

a.  while  VU  is  not  empty 

choose  a  random  vertex  a  from  VC 
choose  a  random  vertex  b  from  VU 
record  an  edge  from  a  to  b 
move  b  from  VU  to  VC 

5.  while  Number  edges  <  Number  desired 
choose  random  vertex  a 
choose  random  vertex  b 
While  edge  (a*b)  exists  or  i f  a  -  b 

increment  a  until  a  -  n  then  increment  b  ana  set  a  -  l 
if  a  *  n  and  b  =  n  then  set  a  -  1  and  b  -  1 
record  edge  (a*b) 

End 
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At  a  later  point  in  the  research*  it  was  decided  that 
for  truly  valid  test  results*  the  graph  generator  must  gen¬ 
erate  purely  random  graphs  with  no  requirement  for  connec¬ 
tivity.  It  was  discovered  that  the  connectivity  requirement 
changed  the  character i st i cs  of  the  sparse  graphs  and  non¬ 
represent  at i ve  graphs  were  being  used.  Therefore*  a  much 
simplified  generator  was  used  which  ignored  steps  1  through 
4  of  GENGRAPH  and  merely  added  random  edges  until  the  re¬ 
quired  number  was  present.  Another  algorithm  called  DIVIDE 
(to  be  described  in  a  later  subsection)  was  then  called  to 
divide  the  graph  into  its  connected  components. 

The  random  number  generator  used  is  identical  to  Grogono 
[Ref.  9J  with  a  modification  to  return  integers  in  any  range 
desired.  The  listing  for  both  GENGRAPH  and  the  Random 
number  generator  is  enclosed  in  appendix  A. 


8.  GRAPH  DIVIDING 


It  has  already  been  shown  that  an  unconnected  graph  nay 
be  colored  by  dividing  it  into  its  connected  components  and 
separately  coloring  each  component.  A  further  simplifies/ 
t i on  can  be  made*  however#  if  even  the  connected  components 
are  divided  at  their  articulation  points#  if  any  exist.  Re¬ 
call  that  an  articulation  point  is  a  vertex  which  is  in 
every  path  from  some  vertex  x  to  some  other  vertex  y.  If  the 
articulation  point  is  divided  into  2  vertices  such  that  the 
oath  from  x  to  y  is  cut#  then  the  graph  is  divided  into 
disconnected  subgraohs.  For  illustration  purposes  consider 


The  articulation  points  are  vertices  d  and  f 


The 


A  Graph  Divided  at  its  articulation  points. 

Figure  o 

In  co'nmunicat  ions  and  t ransport at i on  problems  the  iden¬ 
tification  of  articulation  ooints  is  very  important  because 
these  are  the  'choke*  points?  or  places  where  the  network 
becomes  most  vulnerable  to  failures.  Since  dividing  a  graph 
at  its  articulation  points  results  in  biconnected  components 
which  can  be  colored  separately?  it  is  also  important  to 
discover  the  articulation  points  in  the  colorability  prob¬ 
lem.  The  number  of  colors  needed  will  be  exactly  the  same 
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as  the  number  needed  to  color  tne  component  which  uses  the 
most  colors.  If  the  actual  color  assignments  are  needed, 
then  the  graph  can  be  ’pieced'  back  together  one  component 
at  a  time.  As  each  component  is  added,  the  color  classes 
are  permuted  so  that  the  articulation  vertex  is  in  the  same 
class  in  both  components  of  the  graph. 

An  algorithm  to  find  the  articulation  points  and  the 
biconnected  components  of  a  graph  utilises  a  depth  first 
search  of  the  vertices  in  the  graph  and  maxes  use  of  the 
fact  that  v  is  an  articulation  point  if  and  only  if  there 
are  two  vertices  x  and  y  such  that  every  oath  from  x  to  y 
includes  v. 

The  idea  of  the  depth  first  search  is  to  visit  all  the 
vertices  of  the  graph  bv  starting  at  some  arbitrary  start 
vertex  s  and  recursively  visiting  every  other  vertex.  Given 
that  the  search  is  at  some  vertex  a.  it  follows  an  edge 
(a.b)  to  vertex  b.  If  b  has  already  been  visited.  the 
search  returns  to  a  and  tries  another  edge.  If  b  has  not 
been  visited,  then  the  same  method  is  applied  recursively  to 
b.  At  the  ooint  where  all  the  edges  incident  on  a  have  been 
examined,  the  search  backs  up  along  the  edge  that  took  it  to 
a.  The  search  terminates  at  the  point  where  it  tries  to 
back  uo  from  the  start  vertex.  If  the  graph  was  not  con¬ 
nected#  then  some  new  start  vertex  is  chosen  from  the  un- 
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visited  ones  end  the  procedure  repeated  until  all  vertices 
have  been  visited.  The  edges  which  take  the  search  from  one 
vertex  to  a  vertex  which  has  previously  been  visited  are 
called  'back  edges.'  Edges  which  go  to  an  unvisited  vertex 
are  called  'forward  edges.'  These  definitions  will  be  need¬ 
ed  for  the  algorithm  description. 

Figure  7  shows  a  graoh  with  5  vertices  and  7  edges.  If 
the  depth  first  search  is  started  at  vertex  a  then  the  solid 
arrows  represent  a  possible  set  of  forward  edges  and  the 
dashed  arrows  are  the  back  edges.  The  labels  on  each  edge 
indicate  the  order  in  which  they  are  traversed. 


A  Oeoth  first  search  of  a  Graph. 

Figure  7 

The  basic  idea  of  using  a  depth  first  search  to  find  the 
articulation  ooints  of  a  graph  can  be  seen  by  studying  fig. 
8.  This  schematic  reoresents  a  graoh  with  3  piconnected 
components  labeled  A,  8,  and  C  which  are  joined  at  articula¬ 
tion  vertices  v  and  w.  S  is  an  arbitrary  start  vertex  in  A. 
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A  Schematic  of  a  graph  with  3  biconnected  components  and 
2  articulation  points. 

Figure  8 

If  the  depth  first  search  is  started  at  vertex  S  in  A * 
it  will  eventually  enter  into  C  through  vertex  v.  By  the 
depth  first  nature  of  the  search*  all  the  eoaes  in  C  will  be 
traversed  before  the  search  backs  up  through  v.  If  the 
edges  are  placed  on  a  stack  as  they  are  traversed*  then  when 
the  search  gets  back  to  v*  all  the  edges  in  C  will  be  on  too 
of  the  stack. 

In  order  to  recognize  the  articulation  points*  note  that 
no  vertex  in  C  will  have  a  back  edge  to  any  vertex  which  was 
visited  prior  to  v.  So*  if  the  vertices  are  numbered  as 
they  are  visited  and  each  vertex  is  tagged  with  an  jnae« 
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equal  to  the  smallest  number  of  any  vertex  which  can  be 
reached  throuqh  any  number  of  forward  edges  and  one  back 
edge/  then  the  articulation  points  can  be  found.  When  the 
search  backs  up  from  b  to  a  along  edge  (a#b )/  if  the  index 
tagged  to  b  is  not  less  than  the  number  assigned  a /  then  a 
is  an  articulation  point(or  the  root  of  the  tree).  In  other 
words/  no  vertex  below  a  in  the  tree  has  a  back  edge  to  a 
vertex  visited  before  a.  The  following  algorithm  by 
BaaselRef.  61  was  used  to  find  both  the  connected  components 
and  the  biconnected  components  of  a  graph: 

Algorithm  DIVIDE 


3 V  a  stack  of  vertices 

SE  a  stack  of  edqes 

top  a  top  element  of  SV 

Number  a  array  to  number  the  vertices  in  the  order  visited 
Back  a  array  to  record  the  last  vertex  reachable 
through  a  back  edge. 

Algorithm: 

1.  initialize  Number(i)  =  0/ 

2.  choose  arbitrary  vertex  S 

3.  Number ( 3 )  =  1 

а.  Num  -2  //  next  number 

5.  stack  S  on  3V  (top  =  3) 

б.  while  there  exists  an  unprocessed  edge  from  top  to  v  do 

stack  edge  (too/v)  on  SE 

if  number(v)  >  0  then  back(too)  -  m i n ( bac k ( t op ) / bac k ( v )  ) 

else  do 

number(b)  a  num 
num*e 

back(v)  a  number(v) 
stack  v  on  Sv 
end 

end 
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7.  if  there  is  more  than  one  vertex  on  SV 
let  v  =  second  from  too 
if  back(too)  >s  Number(v)  then 

let  edge  (a*b)  be  the  top  edge  on  SE 

while  Number(a)  >=  Number(v)  and  Number(b)  >= 
Numbe  r ( v  ) 

oop  S£ 

else  Back(v)  =  m i n (Bac k ( v ) * Bac k ( t oo ) ) 
poo  SV 
goto  6  END 

The  code  used  to  impliment  this  algorithm  is  enclosed  in 
append! x  A . 

C.  GRAPH  COLORING  ALGORITHM 

The  essential  idea  of  the  basic  coloring  algorithm  is  to 
build  the  search  tree  until  the  optimal  coloring  is  found. 
Each  node  in  the  tree  represents  the  placement  of  a  vertex 
into  a  color  class  and  the  oath  from  the  root  to  the  node 
reoresents  a  oartial  coloring  of  the  graph.  The  cost  of  a 
node  is  the  number  of  different  colors  used  for  that  partial 
coloring.  As  the  nodes  of  the  tree  are  generated*  by  itera¬ 
tively  expanding  from  the  start  node*  they  are  placed  in  a 
priority  queue  such  that  the  first  node  in  the  queue  is  the 
node  which  reoresents  the  least  number  of  colors  used  and  is 
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deeper  in  the  tree  than  any  other  node  with  the  same  number 
of  colors.  A  node  is  expanded  by  generating  all  of  its  pos¬ 
sible  children  and  the  node  to  be  expanded  next  is  always 
the  first  node  in  the  queue.  The  stopping  point  for  the  al¬ 
gorithm  is  thus  the  point  at  which  the  highest  priority  node 
is  a  full  coloring  of  the  graph.  At  that  point  no  other 
node  in  the  tree  could  possibly  be  expanded  in  any  way  which 
results  in  a  coloring  using  fewer  colors.  The  chromatic 
number  of  the  graph  is  then  known  to  be  the  number  of  colors 
used  in  that  coloring.  Note  that  if  all  possible  colorings 
are  desired  which  use  the  optimal  number  of  colors#  then  the 
algorithm  need  only  save  the  chromatic  number  and  to  contin¬ 
ue  expanding  nodes  until  the  cost  of  the  highest  priority 
node  is  greater  than  the  chromatic  number.  At  each  point 
where  a  node  is  generated  which  is  a  full  coloring,  it  can 
be  printed  out  and  all  optimal  coloring  assignments  will 
thus  be  generated. 

There  are  2  primary  data  structures  used  in  this  algo¬ 
rithm.  These  are  the  TREE  and  QUEUE  data  structures.  The 
elements  of  TREE  are  trees  made  up  of  nodes  and  the  links 
between  them  which  organize  the  nodes  into  a  tree  structure. 
A  node  N  represents  a  color  assignment  to  one  vertex  of  the 
graph.  The  path  from  N  to  the  root  of  the  tree  represents  a 


partial  coloring  of  the  graoh  and  the  cost  assigned  to  node 


N  is  the  number  of  colors  used  in  that  partial  coloring. 
Each  node  has  a  pointer  to  its  oarent  in  the  tree  in  order 
to  maintain  the  links  between  nodes.  There  are  2  operations 
defined  for  TREE.  The  operation  CREATE  generates  a  tree 
with  only  one  node  called  the  root.  The  operation  EXPAND 
takes  as  input  a  tree  T  and  a  particular  leaf  node  N  of  T 
and  returns  a  tree  T  which  is  the  original  tree  T  plus  all 
the  possible  children  of  node  N. 

The  nodes  of  the  tree  are  implimented  using  an  array  of 
records  where  each  record  contains  the  following  5  fields: 

1)  the  depth,  or  level  in  the  tree  at  which  it  is  located, 

2)  the  cost  assigned  to  the  node,  3)  the  vertex  number  which 
was  colored  to  create  the  node,  4)  the  color  assigned  to 
that  vertex,  and  5)  a  pointer  to  the  node's  parent  in  the 
tree.  The  operations  on  the  TREE  are  implimented  in  the  al¬ 
gorithm  CQLQRALL  and  EXPAND  as  follows: 


COLORALL  (Color  all  vertices  of  the  graph) 

1 .  Get  a  node  N 

2.  Get  number  of  first  vertex  -  v  . 

3.  Let  N  be  the  root  node 

N. vertex  =  v 

N.cost  -  l 

N. depth  =  1 

N.eol or  =1 

N. parent  Pointer  a  NULL 

u.  while  N. depth  <  n  do 

expand  N  (see  algorithm  EXPAND  below) 

N  a  first  node  in  priority  queue 
5.  Print  out  3 t a t i s t i c s , co 1 o r i ng, and  chromatic  number. 
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The  algorithm  to  expand  a  node  is  as  follows: 


EXPAND(N)  (N  is  the  node  to  be  expanded) 

1.  Choose  a  vertex  v  to  be  added  to  the  partial  coloring, 
(described  in  subsection  0.) 

2.  let  C  -  { 1 , 2,  .  .  .  » N .cost >  be  the  set  of  colors  used 

in  the  path  to  node  N. 

3.  Trace  the  parent  pointers  from  N  to  the  root.  For 
each  node  M  in  the  oath  to  the  root,  if  an  edge  exists 
between  vertex  v  and  M. vertex  then  remove  color  color 

from  C . 

4.  For  each  color  c^  left  in  C  create  a  new  node  N.  and 

put  it  in  the  priority  queue.  1 

N • .vertex  =  v 

N. .cost  s  N.cost 

N.  .depth  =  N. depth  +  1 

N I . color  =  c • 
n!  .oarpt  r  =  N 

5.  Add  one  new  color  and  create  one  more  node  N.. 

N . .vertex  =  v  J 

NJ. .cost  =  N.cost  *1 
N^. depth  =  N. depth  >1 

.col  or  =  N.cost  +  1 
N  .  .oarpt  r  =  N 
END  J 


In  terms  of  graph  coloring,  step  4  creates  a  child 
node  for  each  of  the  wavs  in  which  the  new  vertex  can  be  ad” 
ded  to  one  of  the  existing  color  classes.  Step  5  generates 
a  node  which  represents  putting  the  new  vertex  in  a  new 
color  class  by  itself. 

The  QUEUE  data  structure  is  used  to  organize  the  leaf 
nodes  of  a  tree  according  to  their  priority.  There  are  2 
components  which  determine  the  priority  of  a  node.  The 
first  is,  of  course,  its  cost  as  defined  earlier.  The 
second  component  is  used  to  determine  the  highest  priority 
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between  2  nodes  with  the  same  cost.  nhen  more  than  one  node 
is  available  which  have  the  same  cost  then  the  node  which  is 
deeoest  in  the  tree  has  a  higher  priority  since  it  is  closer 
to  a  complete  coloring.  There  are  two  operations  defined 
for  the  QUEUE  called  ADDNODE  and  REMOvETOP.  ADDNQDE  takes  a 
queue  Q  and  a  node  N  and  returns  a  queue  Q  with  N  in  its 
proper  priority  location.  REMOVETOP  takes  a  queue  Q  and  re¬ 
turns  the  node  N  with  highest  priority. 

The  priority  queue  is  implimented  as  a  heap*  as 
described  in  ref.  10.  The  nodes  are  organized  into  a  full 
binary  tree  such  that  the  root  is  the  node  with  highest 
priority  and  every  node  in  the  tree  has  a  higher  priority 
than  either  of  its  children.  When  adding  a  node  to  the 
t  ree,  it  is  put  in  the  first  available  location  and  then 
'sifted  up'  until  its  parent  has  a  higher  priority.  When 
the  root  node  is  removed*  the  last  node  is  moved  to  the  root 
position,  then  'sifted  down'  until  both  its  children  have  a 
lower  priority.  The  binary  tree  is  maintained  by  an  array 
of  pointers  qO  where  q ( 1 )  Points  to  the  node  at  the  root  of 
the  binary  tree  and  the  children  of  q(i)  are  located  at 
q(2i )  and  o'2i  ♦  1 )  . 

In  any  tree  search  algorithm  it  is  important  to  ensure 
that  isomorphic  paths  are  not  generated.  For  the  problem  of 
qraoh  co 1  or ab i 1 i t y ,  two  colorings  are  isomorphic  if  the  ver- 
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tices  of  the  graph  are  partitioned  into  color  classes  in  the 
same  way  in  both  colorings.  Therefore#  if  one  coloring  uses 
k  color  classes#  then  there  are  k!  possible  isomorohs  of 
that  coloring.  In  order  to  prevent  the  tree  from  containing 
all  of  these  unneeded  paths#  one  must  either  check  for  pos¬ 
sible  isomorohs  with  the  addition  of  each  new  node  or 
prevent  isomorphic  path  generation  from  occuring  by  the 
design  of  the  algorithm. 

There  are  two  essential  ingredients  of  algorithms 
COLORALL  and  EXPAND  which  ensure  no  isomorphs:  1)  When  a 
node  N  is  expanded#  at  most  1  new  color  is  added  to  the 
available  color  set#  and  2}  In  expanding  a  node  N#  only  1 
vertex  of  the  uncolored  vertices  can  be  added  to  the  partial 
coloring  and  that  vertex  must  be  the  same  for  all  the  chil¬ 
dren  of  N. 

THEOREM:  Any  tree  search  graph  coloring  algorithm  with  the 
above  2  properties  will  not  generate  isomorphic  colorings. 

Proof:  Assume  the  algorithm  generates  a  tree  with  2  iso¬ 
morphic  paths  P  and  Q.  Then,  by  the  definition  of  isomorph¬ 
ic  colorings#  for  each  color  class  in  P  there  exists  a  color 
class  in  Q  with  exactly  the  same  elements.  Let  vertex  v  be 
the  first  vertex  colored  in  P  and  Q  such  that  the  color  of  v 
in  P  is  different  from  the  color  of  v  in  Q.  Call  the  class 
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p 

containing  v  in  path  P*  C  j  and  the  class  containing  v  in  Q* 
q 

Cj  •  If  the  subscript  indicates  the  color*  then  i  t  j. 

Since  only  one  new  color  could  be  added  when  v  was  colored* 

P  q 

then  at  least  one  of  Cj  or  Cj  contains  a  vertex  w  which  was 
colored  prior  to  v.  without  loss  of  generality*  let  that 
class  be  C?.  In  other  words*  v  and  w  are  in  the  same  color 
class  in  path  P.  But*  since  v  was  colored  differently  in 
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path  Q  and  w  is  the  same  in  both  P  and  Q*  there  can  be  no 
color  class  in  Q  containing  both  v  and  w.  Therefore*  P  and 
Q  are  not  isomorphic  colorings.  Figure  9  illustrates  this 
proof . 


P  Q 

Two  Paths  P  and  Q  which  are  not  isomorphic 

F i gure  9 


D.  THE  VERTEX  ORDERING  FUNCTIONS 

The  purpose  of  this  research  was  to  investigate  dif¬ 
ferent  branching  functions  in  a  branch  and  bound  approach  to 
solving  the  chromatic  number  problem.  For  the  particular 
branch  and  oouno  algorithm  used#  the  branching  function  is 
the  order  in  which  the  vertices  are  selected  to  be  added  to 
the  partial  coloring  at  each  new  level  of  the  tree.  In  this 
research#  five  different  methods  of  ordering  the  vertices 
were  investigated.  Four  of  the  methods  were  static  order¬ 
ings  in  which  the  vertices  were  preoraered  and  then  the 
graph  coloring  algorithm  was  called  to  search  the  entire 
tree.  In  every  oath  through  the  tree  in  these  methods  the 
order  of  choosing  the  next  vertex  was  constant.  In  the 
fifth  method#  the  vertex  chosen  to  be  colored  next  was 
dependent  on  the  particular  color  assignments  made  prior  to 
the  node  to  be  expanded.  This  method  of  dynamic  ordering 
results  in  a  larger  time  expenditure  at  each  node  since  the 
order  must  be  recomputed  at  each  expansion  of  the  tree. 

Note  that  the  purpose  of  any  of  these  orderings  is  to 
make  the  search  tree  as  small  as  possible#  thus  requiring 
fewer  nodes  and  fewer  expansions.  If  the  amount  of  work 
done  to  expand  a  node  (create  its  children)  is  constant  for 
each  method#  then  the  resulting  size  of  the  search  tree  is  a 
direct  measure  of  the  benefit  of  any  procedure.  Since  the 


45 


height  of  the  tree  will  always  be  the  number  of  vertices  in 


the  graph#  the  only  way  to  decrease  the  number  of  nodes  is 
to  decrease  the  branching  factor  (number  of  children  of  each 
node).  The  optimal  branching  factor  is#  of  course*  1.  With 
the  algorithm  COLORALL*  this  is  reached  when  each  vertex  to 
be  colored  requires  a  new  color.  It  is  approached  most 
closely  when  the  vertices  are  ordered  such  that  the  next 
vertex  is  always  the  one  with  the  fewest  colors  available  to 
be  colored.  Three  of  the  static  orderinqs  attempt  to  accom- 
olish  this  ordering  through  the  use  of  some  heuristic  func* 
tion  to  order  the  vertices.  The  dynamic  ordering  actually 
determines#  at  each  node#  which  of  the  uncolored  vertices 
has  the  fewest  colors  available. 

Obviously*  any  graph  which  is  completeCa  clique)#  will 
be  colored  with  a  branching  factor  of  one#  no  matter  what 
order  is  used.  The  same  applies  to  a  graph  which  is  just  a 
forest  of  verticesCno  edges).  The  difficulty  lies  in  the 
midrange.  Another  way  to  look  at  it  is  that  coloring  is 
easiest  when  there  are  few  rest r i c t i ons ( edges )  or  very  many 
restrictions  to  the  colors  which  can  be  used.  Each  of  the 
ordering  methods#  therefore#  really  attempts  to  find  the 
vertices  with  the  most  restrictions  and  color  them  first. 
Thus#  the  tree  has  the  smallest  branching  factor  at  the  too 
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and  more  of  the  tree  is  pruned  off  resulting  in  fewer  nodes 


to  search. 

Following  is  a  description  of  each  ordering.  The  list¬ 
ing  for  all  orderings  is  in  the  program  QRDERNODE  in  appen¬ 
dix  A.  The  case  numbers  in  the  switch  statement  correspond 
to  the  numbers  shown  here. 

1 .  Random 

This  method  of  selecting  vertices  was  used  as  a 
baseline  against  which  to  compare  any  benefit  of  the  other 
methods.  Inasmuch  as  the  graph  was  generated  randomly,  no 
real  work  needed  to  be  done.  The  vertices  were  merely 
Chosen  in  the  order  of  increasing  vertex  number. 

2.  Oegree 

This  method,  used  in  the  welsh-Powell  approximation 
al gor i t hm  [Ref .  Ill,  orders  the  vertices  by  decreasing  values 
of  the  degree  of  the  vertices.  Recall  that  the  degree  of  a 
vertex  is  the  number  of  edges  incident  upon  it  and  is  easily 
calculated  by  summing  the  elements  of  the  rows  of  the  adja¬ 
cency  matrix.  The  idea  behind  this  method  is  that  the  ver¬ 
tex  with  the  highest  degree  has  more  restrictions  on  it  when 
choosing  a  color  so  it  should  be  colored  first.  That  color 
used  is  then  removed  from  the  colors  available  for  every 
vertex  connected  to  it. 
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.  Ove 

OVEC  3  is  the  name  given  by  this  author  to  a  method 
devised  by  M.R.wi II i ams  IRef .  41.  In  the  prior  method  of 
ordering  the  vertices  by  degree*  the  underlying  premise  is 
that  if  a  vertex  i  is  connected  to  p  other  vertices*  then  it 
will  be  more  difficult  to  color  than  vertex  j  which  is  only 
connected  to  g  vertices  (p  <  q).  What  that  method  fails  to 
consider  is  the  degree*  or  difficulty*  of  the  vertices  to 
which  i  and  j  are  connected.  For  example*  if  vertex  j  is  in 
a  clique  of  size  4  and  i  is  the  center  of  a  star  with  5 
points*  as  shown  in  fig.  10*  then  certainly  j  should  be 
colored  before  i*  as  well  as  all  the  other  vertices  in  the 
C 1 i que. 


Example  of  Shortcoming  of  Degree  method 
Figure  10 


Williams  recognized  that  the  vector  of  degrees  of  the  ver¬ 
tices  was  the  first  step  towards  the  development  of  the  dom¬ 
inant  eigenvector  of  the  adjacency  matrix.  Recall  that 


•  z  ... 


J 


1  J 


48 


In  the  standard  technigue* 


d^  is  the  eigenvector 


cor respondi ng  to  the  largest  eigenvalue  of  matrix  A  and  is 
found  by  iteration  using: 


n  +  1 


-  z 

j 


aij* 


If  d.j  is  the  vector  of  degrees*  then  this  scheme 
does  take  into  account  the  degrees  of  the  vertices  to  which 
a  vertex  is  connected.  And*  in  fact*  the  number  of  itera¬ 
tions  used  determines  how  far  away  from  vertex  i  to  iook  in 
calculating  the  relative  difficulty  of  coloring  vertex  i. 

Williams  also  determined  that  the  relative  order  of 
after  a  few  iterations.  This  fact  was  borne  out  by  results 
of  this  research.  Although  improvement  was  seen  out  to  2  or 
3  iterations*  little  or  no  changes  in  the  size  of  the  search 
tree  was  found  after  3  iterations.  The  only  effect  was  the 
increased  computation  to  calculate  the  vectors.  For  this 
reason*  3  iterations  was  chosen  as  a  good  numoer  to  use  to 
compare  this  method  against  the  others. 

4 .  Ad  j  acenc v 

The  rationale  behind  this  ordering  method  is  that 
the  vertex  which  is  adjacent  to  the  most  previously  colored 
vertices  will  probably  have  the  fewest  colors  available. 


i 
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This  is  not  necessarily  true*  however,  as  a  vertex  nay  be 
adjacent  to  many  vertices  which  are  all  in  the  same  color 
class/  thus  reducing  its  available  color  space  by  only  one 
color. 

The  algorithm  to  determine  which  vertex  is  next/ 
considers  the  rows  in  the  adjacency  matrix  for  all  the  ore* 
viously  colored  vertices/  and  adds  the  columns.  The  highest 
column  total  then  indicates  which  vertex  is  to  be  colored 
next.  Previously  colored  vertices  can  be  eliminated  dv  sub¬ 
tracting  an  arbitrarily  high  number  from  the  running  column 
total/  for  each  vertex,  when  it  is  chosen.  In  the  case  of 
ties,  the  vertex  of  highest  degree  is  chosen. 

5.  Least  colors  available  C L C A ) 

In  all  the  previous  static  orderings/  an  assumption 
was  made  that  the  best  vertex  to  color  next/  at  any  point  in 
the  tree/  was  independent  of  the  specific  color  asignments 
made  up  to  that  point.  he  were  simply  trying  to  find  the 
vertex  that  was  'most  probably'  a  better  choice  than  the 
others.  In  this  dynamic  ordering,  when  expanding  each  node/ 
we  actually  choose  the  vertex  which  will  result  in  the 
fewest  children  being  created/  thus  minimizing  the  branching 
factor.  Since  a  child  node  is  created  for  each  available 
color  which  the  next  vertex  can  be  colored/  we  must  consider 
all  uncolored  vertices  in  the  oath  to  the  node  being  expand- 
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ed  and  determine  which  vertex  has  the  fewest  colors  avail 


aole.  To  do  this  note  that  ideally  we  would  like  to  choose 
a  vertex  which  is  adjacent  to  at  least  one  vertex  in  each 
color  class  used  so  far.  Thus  a  new  class  would  have  to  be 
created  and  there  would  be  only  one  child  node  created. 
Short  of  this  ideal  situation,  we  would  like  to  choose  the 
vertex  which  is  adjacent  to  at  least  one  vertex  in  most  of 
the  color  classes  so  far  created. 

This  is  accomolished  by  creating  a  new  matrix  CC 
simitar  to  the  adjacency  matrix  such  that  there  is  a  row  in 
CC  for  each  color  class  generated  and  a  column  for  each  ver¬ 
tex.  Then  we  define 


CCij  * 


j  1  if  (v,i)  6  E  for  some  vertex  v  in  color  class  i 
f  0  otherwise. 


i.e.,  ccij  s  1  iff  vertex  j  cannot  be  colored  color  i. 

The  matrix  CC  is  easily  created  from  the  adjacency 

matrix  A  as  fol lows:  Let  cc  be  row  i  in  LC#  and  a  be  row  i 

i  i 

in  A,  both  represented  as  bit  strings. 

initialise  all  cc.  =  0 

for  i  from  1  to  NU^-COLOR-CL ASSES 


for 

each  vertex  j 

i  n 

color  class  i 

II 

U 

u 

CCj  V 

ai 

The 

ne*t  step  is 

i 

t  0 

generate  a  vector 

i 

X  cc 

i  f  i 

is  uncolored 

SUM.s  j 

i 

’  j 

J 

1  f 

j  i  s 

colored 

SI 


Then  the  value  of  j  for  which  SUM  is  a  maximum  is 
the  vertex  which  has  the  fewest  color  classes  available  to 
which  it  can  be  aaded.  Again*  for  consistency*  ties  are 
broken  by  choosing  the  vertex  of  highest  degree. 

E.  CONTROLLING  ALGORITHM 

Inasmuch  as  the  ourpose  of  this  research  was  to  experi¬ 
ment  with  various  orderings  of  the  vertices*  the  flow  of 
control  was  fairly  st rai qht f orward.  After  only  a  few  trial 
cases  it  was  discovered  that  the  primary  limitation  to  the 
size  of  the  graph  which  could  be  colored  was  computer 
storage  soace*  not  time.  As  the  graphs  got  larger*  more 
nodes  were  required  in  the  search  tree  until  the  availaDle 
memory  was  used  up  at  about  4000  nodes.  Although  unsigned 
variables  were  used  in  the  node  structures  and  space  was 
conserved  wherever  possible*  this  limitation  could  not  be 
s i qn i f i cant  1 v  changed.  For  even  the  best  orderings*  this 
meant  graphs  of  about  30  vertices  or  less.  For  empirical 
data  then,  graphs  were  generated  with  from  10  to  30  vertices 
and  edges  from  20X  to  SOX  of  the  maximum  possible.  For  each 
graph  size,  100  graphs  were  generated  (some  of  which  were 
probably  identical  and  many  which  were  at  least  isomorphic) 
and  each  ordering  method  was  used  to  color  the  graph.  The 
graph  was  then  divided  at  its  articulation  points  and  each 
biconnected  subgraph  was  colored  using  each  ordering  func- 
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tion.  Data  about  the  size  and  shape  of  the  search  tree  was 
collected  and  tabulated  for  all  the  above  samples 
( append ix  8 ) . 

Note  that  100  trials  of  each  graph  size  is  a  relatively 
small  sample#  but  was  a  necessary  restriction  to  keep  the 
entire  program  execution  time  to  within  a  few  hours. 

The  following  algorithm  shows  the  setup  used  for  obtain¬ 
ing  tne  experimental  data. 


al gori thm  MAIN 


initial i z  e 

For  n  from  minsize  to  maxsize  (n  is  number  of  vertices) 
begi  n 

compute  maxedges  =  (n)(n-l)/2 

For  fullness  from  30%  to  70%  of  maxedges 

begin 

compute  number  of  edges  e  =  fullness  *  maxedges 
For  trials  from  1  to  100 
begi  n 

generate  random  graph  (n  vertices#  e  edges) 

For  each  ordering  method 
begi  n 

color  the  graph 

record  data  (expansions#  nodes#  branching  factor) 
divide  the  graph 

For  each  biconnected  component  of  the  graph 
begin 

color  the  component 
record  component  data 

end 

sum  the  component  data 

end 

end 

compute  averages  of  all  data  collected 
print  out  data 

end 

END 
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The  function  MAIN  in  appendix  A  is  the  controlling  algo* 
rithm  for  the  entire  orocess. 

F.  GENERAL  COMMENTS  ON  ALGORITHM  DEVELOPMENT  AND  THE 

PROGRAMMING  LANGUAGE  'C* 

In  all  the  programming  done  during  this  research,  ef¬ 
forts  were  made  to  save  time  and  storage  space  whenever 
feasible.  The  first  attempt  to  design  a  coloring  algoritnm 
used  a  backtracking  method  which  is  common  for  tree 
searches.  However,  the  size  of  the  search  tree  was  so  large 
that  the  storage  required  to  stack  the  variables  used  in 
each  recursive  call  to  the  backtracking  function  was  exorbi¬ 
tant.  For  that  reason,  an  iterative  scheme  was  chosen  and 
i mo  1 i mented . 

Further  savings  were  made  by  using  bit  strinq3  and  1  bit 
variables  whenever  possible.  Since  the  maximum  qraph  size 
was  to  be  less  than  3 2  vertices,  the  long  inteqer  variables 
of  32  bits  in  *C'  were  ideally  suited  for  the  adjacency  ma¬ 
trix  and  many  other  needed  variables.  Also  'C'  has  a  good 
assortment  of  bit  manioulation  operators  including  multiple 
shifts,  or,  and,  and  exclusive  or.  Using  these  led  to  time 
savings  in  many  cases. 

Space  was  conserved  by  using  the  unsigned  variable  type 
in  arrays  and  records.  With  this  convention,  the  length  of 
the  storage  cells  reserved  were  only  as  long  as  needed.  For 
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example*  the  vertex  numbers  from  0  to  31  could  be  stored  in 
only  5  bits.  Another  oract i ce<  not  necessarily  considered 
good  orogramminq  technique  was  the  fairly  extensive  use  of 
global  variables.  Though  risky*  «t  did  save  the  time  re¬ 
quired  to  pass  parameters  and  the  space  to  store  local  vari¬ 
ables. 

Overall*  this  author  feels  that  *  C  *  is  a  very  appropri¬ 
ate  language  for  research  of  this  type. 


IV.  RESULTS 


A .  THE  PARAMETERS 

In  attempting  to  determine  the  relative  'goodness'  of 
several  computer  programs#  one  must  first  find  some  measur- 
aole  parameters  to  be  minimized  or  maximized  to  indicate  the 
best  orogram.  For  approximation  algorithms#  one  is  very 
often  not  Interested  in  how  long  a  program  takes  as  long  as 
it  is  some  linear  function  of  the  Droblem  size  or  possioly 
even  O(nlgn).  What  one  measures  is  the  closeness  of  the  ap¬ 
proximation  to  the  actual  optimal  solution. 

For  this  research#  the  solution  to  the  graph  coloring 
problem  was  to  be  the  exact  chromatic  number  in  every  case. 
The  coloring  algorithm  Mas  to  find  the  chromatic  number  of  a 
given  graph  and  the  'goodness'  of  the  vertex  ordering 
methods  was  to  be  determined  by  how  fast  the  algorithm  ran 
and  how  much  space  was  used.  The  specific  questions  to  be 
answered  were:  1)  Is  there  a  best-time  ordering  for  all 
graphs?#  2}  Is  there  a  best-soace  ordering  for  all  graphs?# 
3)  Are  they  the  same?#  What  type#  or  size#  of  graphs 
favor  which  ordering  if  there  is  a  difference?#  and  5)  Is 
one  ordering  better  on  the  average#  but  another  better  for 
worst-case  situations. 

The  results  of  the  many  experiments  which  relate  to 
these  questions  is  tabulated  in  appendix  B.  The  discussion 
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here  is  of  a  general  nature  and  will  summarize  the  results. 
There  were  basically  two  separate  categories  of  measurements 
made  -  those  related  to  time  and  those  related  to  space. 
Though  they  are  closely  related  in  this  problem,  the  parame¬ 
ters  measured  are  distinct  and  will  be  discussed  separately. 

1 .  T  j  me 

In  determining  the  difference  in  time-complexity  of 
the  5  coloring  methods,  there  are  2  modules  of  the  program 
to  consider  -  ORDERNODE  and  CQLORALLCof  which  EXPAND  is  a 
part).  These  functions  are  the  only  two  whose  time- 
complexity  is  a  function  of  the  vertex  ordering  as  well  as 
the  graph  size. 

For  the  4  static  orderings  the  time-complexity  is 
easy  to  determine.  For  consistency,  each  ordering  merely 
filled  an  array  called  next,  where  next(i)  was  the  vertex  to 
follow  i  in  the  ordering.  Random  is  therefore  clearly  0(n) 
since  the  vertices  are  ordered  bv  vertex  number.  One 

traversal  through  the  array  next  is  all  that  is  needed  to 

2 

fill  it.  DEGREE  is  0(n  )  since  the  adjacency  matrix  has  to 
he  traversed  to  find  the  degree  of  each  vertex,  then  a  sim- 
ole  sortinq  of  the  vertices  is  O(nlgn)  so  the  overall  com- 
plexity  is  0 ( n  ) •  DVEC  3  is  0  ( n  )  if  the  number  of  itera¬ 
tions  is  kept  const  ant ( such  as  3  in  this  case).  It  would  be 
conceivable  to  iterate  n  times  for  the  best  possible  order- 
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ing  in  which  case  DVEC  n  would  be  0(n  ).  ADJACENCY  is  also 

2 

0(n  )  because  a  running  total  was  kept  of  the  relative  oosi- 
t i on  of  the  unselected  vertices  as  each  new  one  was  deter¬ 
mined.  Thus  the  adjacency  matrix  was  addressed  in  only  one 
row  for  each  new  vertex.  For  the  dynamic  ordering  LCA.  no 
work  was  done  in  the  function  ORDERNODE  but  rather  in  the 
function  PICKNODE  which  was  called  every  time  a  node  was  to 
be  exoanded.  PICKNODE  itself  has  a  time  comolexity  of  n^ 
since  it  takes  n4  operations  to  create  the  n  x  n  matrix  of 
color  classes  CC  and  each  element  in  CC  reauires  0(n)  to 
compute.  But  the  bigger  concern  is  that  the  time -complexity 
of  the  branching  function  EXPAND  is  n^  times  as  great  for 
LCA  than  for  the  static  orderings. 

For  the  static  orderings  the  time-complexity  of  gen¬ 
erating  the  search  tree  has  2  components  -  the  time  to  ex¬ 
pand  a  node  and  the  number  of  expansions  needed  to  find  a 
solution.  For  LCA/  a  third  component/  the  time  to  find  the 
next  vertex,  must  also  be  considered. 

Once  the  next  vertex  is  found/  the  time  reouired  by 
the  function  EXPAND  is  0(n).  The  only  iterative  part  of 
that  function  is  following  the  oath  to  the  root  and  then 
generating  a  child  node  for  each  possible  coloring.  Since 
the  depth  of  the  tree  can  be  up  to  n-1  and  the  number  of 
children  can  be  n.  this  function  is  0(n). 
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So  far»  all  the  time  complexities  have  been  polyno 


mial.  If  this  were  also  true  of  the  number  of  expansions* 
then  one  could  claim  that  the  overall  algorithm  would  run  in 
doI ynomi al -t i me  since  the  product  of  any  number  of  polynomi¬ 
als  is  also  a  polynomial.  Unfortunately*  that  is  not  the 
case  for  worst-case  situations.  Consider  for  a  moment  the 
sire  of  the  full  search  tree  generated  by  this  algorithm  if 
no  pruninq  is  allowed.  Refer  to  figure  11  which  shows  the 
first  four  levels  of  the  search  tree.  The  notation  x*y  in 
each  node  of  the  tree  indicates  that  the  node  represents  an 
assignment  of  color  y  to  vertex  x.  Note  that  at  the  root 
there  is  t  node  followed  by  2  nodes  at  level  2.  Then  at 
level  three  there  are  5  nodes*  2  children  of  one  in  level  2 
and  3  for  the  other.  To  determine  the  number  of  nodes  at 
level  k *  recall  that  one  color  class  is  added  for  each  level 
in  the  tree.  Also*  a  node  at  level  k  represents  one  of  tne 
ways  of  partitioning  the  vertices  into  k  or  fewer  color 
classes.  Therefore*  the  number  of  nodes  at  level  k  is  the 
number  of  ways  to  partition  a  set  of  n  elements  into  k  or 
fewer  classes. 

The  Stirling  numbers  S(n*m)  represent  the  numDer 
of  ways  to  partition  n  elements  into  m  distinct 
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THE  FIRST  FOUR  LEVELS  OF  THE  UNPRUNED  SEARCH  TREE 


classeslRef.  12)  and  are  generated  by  the  following  recur¬ 
rence  relation: 

S(n*l*k)  -  S(n,k-l)  ♦  k  *  S(n,k) 

where  S  ( n  *  n  )  -  1  \ 

/  for  all  n 

and  S ( n / 1 )  -  1  ) 

To  understand  this  relation*  consider  the  following: 
Given  that  n  elements  can  be  partitioned  into  k  classes  in 
S ( n  *  k )  ways*  there  are  two  ways  to  add  the  (ntl)st  element. 
First*  it  can  be  added  to  any  one  of  the  k  classes  for  each 
way  of  partitioning  n  elements.  The  number  of  ways  of  doing 
this  is  given  by  k  *  S(n*k).  Secondly*  the  n  elements  can 
be  Partitioned  into  k-1  classes  and  then  the  (ntl)st  element 
can  be  a  single  element  in  a  new  class.  This  is  the  reason 
for  the  S(n*k-1)  term. 

But  the  number  of  nodes  at  level  k  is  the  sum  of  the 
Stirling  numbers  S(n*m)  for  m  from  1  to  k.  These  numbers 
are  called  the  Bell  numbers  9  after  E . T  .  Be  1 1  [Re f  .  13).  B 

n 

is  the  number  of  partitions  of  a  set  with  n  elements  into 
any  number  of  distinct  classes  up  to  n.  That  is 

Bn  *  £  S(n, k) 

k 

To  see  how  fast  these  numbers  grow*  the  first  8  are: 
1,2*5*15, 52*203*877, aiao 

In  other  words*  at  level  8  in  the  search  tree  there 
would  be  4140  nodes  and  the  total  numoer  of  nodes  in  the 


tree  would  be  the  sum  of  8 ..  for  i  from  1  to  8.  To  see  what 
the  uoper  and  lower  bounds  on  Bn  are#  consider  again  the 
tree  in  fig.  11.  The  fewest  number  of  children  of  each  node 
is  found  in  the  leftmost  path  down  the  tree.  In  this  path# 
every  node  has  two  children.  In  the  rightmost  oath  down  the 
tree#  a  node  at  level  i  has  i  children  and  this  is  the  larg¬ 
est  number  of  children  for  any  node  at  that  level.  Thus  the 
total  number  of  nodes  at  any  level  i  lies  somewhere  between 
the  two  cases#  or 


Since  Bn  is  growing  faster  than  2n  it  is  easy  to  see 
that  as  n  gets  large#  the  possible  size  of  the  tree  is  going 
to  far  outweigh  any  of  the  polynomials  so  far  discussed. 
Obviously#  unless  a  significant  amount  of  pruning  is  done# 
the  time  to  expand  every  node  at  each  increasing  level  in 
the  tree  is  going  to  grow  far  too  rapidly.  Since  the  search 
strategy  was  fixed  for  all  orderings#  the  most  significant 
way  to  compare  the  time  savings  of  each  orderinq  was  to 
count  the  number  of  nodes  which  needed  to  De  expanded. 
Therefore#  the  number  of  calls  to  the  function  EXPAND  was 
used  as  the  primary  measure  of  the  time  complexity  of  the 


various  orderings 


2 .  Space 

There  are  only  two  data  structures  in  these  programs 
whose  sile  is  dependent  upon  the  vertex  ordering:  the  array 
of  nodes  in  the  search  tree  and  the  array  of  pointers  used 
for  the  priority  queue.  But  the  number  of  nodes  in  the 
queue  is  always  less  than  the  total  number  of  nodes  in  the 
tree  since  expanded  nodes  are  removed  from  the  queue.  Since 
the  number  of  expansions  has  been  recorded  for  a  time 
analysis*  only  the  total  number  of  nodes  in  the  search  tree 
at  the  time  an  optimal  solution  is  found  need  be  recorded 
for  a  space  analysis. 

In  order  to  see  the  shape  of  the  search  tree*  anoth¬ 
er  parameter*  the  average  branching  factor*  was  recorded  for 
each  level  in  the  tree.  The  average  branching  factor  is  de¬ 
fined  to  be  the  average  number  of  children  of  each  node  in 
the  search  tree.  A  full  binary  tree  would  thus  have  an 
average  branching  factor  of  2.0.  If  the  branching  factor  is 
recorded  for  each  level  in  the  search  tree*  one  can  see  the 
'shape'  or  the  'narrowness'  of  the  tree.  The  most  desirable 
situation  is  the  one  in  which  the  branching  factor  stays  at* 
or  as  close  as  possible*  to  1.0  in  the  first  several  levels 
of  the  tree.  This  indicates  that  the  maximum  pruning  took 
place  as  early  in  the  search  as  possible*  thus  reducing  the 
search  space  by  more  nodes  for  each  prune.  The  'goodness' 
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of  the  ordering  functions  can  thus  be  seen  by  how  long  the 
branching  factor  stays  at  1.0  and  how  low  it  stays  when  it 
does  start  to  increase. 

B.  TIME  ANALYSIS 

In  analysing  the  time  required  by  the  program  for  each 
ordering  method  there  were  two  cases  of  interest:  average 
time#  and  the  time  required  in  the  worst  case.  Figure  12  is 
a  plot  of  the  average  number  of  expansions  required  for  each 
ordering  method  when  graphs  containing  30%  of  the  maximum 
edges  were  colored. 

For  the  worst  case  graphs  the  natural  logarithm  of  the 
number  of  expansions  was  plotted  since  exponential  growth 
was  expected.  Figure  13  is  a  plot  of  the  worst  case  graphs 
for  30%  edges.  These  plots  are  not  smooth  and  in  fact  are 
quite  ragged  because  each  point  actually  represents  the 
coloring  of  one  particular  graph  in  the  sample  of  100.  Note 
that  the  ordering  methods  parallel  each  other  in  the  dies 
and  hills  of  the  plot*  indicating  that  the  same  graph  was 
usually  the  most  difficult  to  color  for  all  of  the  methods. 

The  significance  of  these  two  plots  is  in  the  total  dom¬ 
inance  of  the  LEAST  COLOR  AVAILABLE  method  of  ordering  the 
vertices.  Though  all  the  other  methods  eventually  'explod¬ 
ed'  even  in  the  average  case*  LCA  stayed  very  close  to 
linear  for  the  full  range  of  graphs  up  to  32  vertices. 


AVERAGE  EXPANSIONS  VS.  GRAPH  SIZE 
Figure  12 


LN  (MAXIMUM  EXPANSIONS) 


Furthermore*  the  worst  case  analysis  shows  the  number  of  ex 


oansions  for  LCA  leveling  off  between  20  and  30  vertices. 
It  would  be  an  interesting  experiment  to  carry  the  graph 
sires  out  past  32  to  see  if  and  when  LCA  also  starts  to  show 
explosive  qrowth. 

Another  result*  though  not  shown  in  graphical  form  is 
that  the  time  required  by  any  of  the  ordering  methods  is 
lowest  for  graphs  with  very  few  edges  as  well  as  for  graphs 
which  are  close  to  being  complete.  For  any  given  value  of 
n*  empirical  evidence  indicated  that  the  most  time  is  re¬ 
quired  when  the  number  of  edges  is  in  the  vicinity  of  50% 
(see  appendix  B).  This  is  an  intuitively  obvious  result  if 
one  considers  what  happens  in  the  search  tree.  fv  i  t  h  very 
few  edges*  very  few  colors  will  be  needed  so  the  number  of 
children  of  a  node  will  be  minimal.  The  search*  in  this 
case*  will  tend  to  proceed  down  the  left  side  of  the  search 
tree.  Most  of  the  right  side  will  be  pruned  away  Because 
the  cost  of  the  nodes  in  the  right  side  of  the  tree  will  be 
greater  than  the  cost  of  the  optimal  solution. 

In  the  case  where  there  are  many  edges*  a  new  color  will 
be  required  at  each  level  and  the  search  will  tend  to  stay 
very  close  to  the  right  side  of  the  tree.  The  left  side 
will  be  pruned  away  since  many  of  the  nodes  on  that  side 
will  represent  illegal  colorings  due  to  edaes. 
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C.  SPACE  ANALYSIS 


Just  as  for  time/  the  space  analysis  is  of  Doth  average 
and  worst  case  situations.  Figure  14  is  a  plot  of  the  aver¬ 
age  number  of  nodes  in  the  search  tree  when  the  first  op¬ 
timal  solution  was  found.  Figure  15  is  a  plot  of  the  natur¬ 
al  log  of  the  maximum  nodes  needed  for  any  of  the  100 
graphs.  Again/  the  number  of  edges  is  30%  of  the  maximum 
possible. 

Not  surprisingly/  the  space  analysis  results  are  very 
similar  to  the  results  found  in  the  time  analysis.  In  every 
case  LC A  is  superior  in  the  sense  that  it  uses  less  space 
and  its  growth  rate  is  much  improved  over  the  other  metnods. 

Since  the  relationship  between  the  number  of  expansions 
and  the  total  number  of  nodes  is  the  branching  factor,  it 
comes  as  no  surprise  that  LCA  exhibits  the  lowest  branching 
factor  at  the  too  of  the  search  tree.  Figure  16  shows  the 
average  branching  factor  as  a  function  of  depth  for  each  of 
the  4  heuristic  orderings.  Note  that  LCA  and  ADJACENCY  have 
the  same  branching  factor  for  3-4  levels  but  then  ADJACENCY 
grows  faster  in  the  midrange.  Note  also  that  both  of  these 
methods  have  a  higher  branching  factor  in  the  middle  of  the 
tree  than  either  DEGREE  or  DVEC  3.  This  fact  shows  how 
powerful  3  or  4  pruninqs  at  the  top  of  the  tree  can  be. 
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AVERAGE  NODES  VS.  GRAPH  SIZE 
Figure  14 


MAXIMUM  NODES  VS.  GRAPH  SIZE 
F i gu  re  15 


AVERAGE  BRANCHING  FACTOR  VS.  DEPTH  IN  TREE 
Figure  16 


0.  ADVANTAGE  OF  DIVIDING  THE  GRAPHS 


In  the  first  stage  of  this  colorinq  process#  the  grapns 
were  colored  just  as  generated  with  no  attempt  to  divide 
them  into  biconnected  components.  In  the  second  staqe#  they 
were  divided  and  then  each  component  was  colored  and  the 
number  of  expansions  required  was  summed  up  for  all  the  com¬ 
ponents.  For  sparse  graphs  where  there  were  many  very  small 
components#  the  results  of  this  experiment  were  quite  sub¬ 
stantial.  However#  as  the  graphs  became  more  dense#  the 
original  graph  was  biconnected  and  no  savings  were  derived 
f  rom  dividing. 

Figure  17  shows  the  size  of  the  largest  biconnected  com¬ 
ponent  as  a  function  of  the  number  of  edges.  In  this  plot# 
the  ratio  of  the  size  of  the  largest  component  to  the  origi¬ 
nal  graph  size  is  plotted  for  4  cases.  Note  that  since  the 
maximum  number  of  edges  grows  as  n^  ,  when  the  number  of 
vertices  in  the  graph  increases#  the  graph  becomes  bicon¬ 
nected  with  a  smaller  percentage  of  edaes.  Thus  it  is  rea¬ 
sonable  that  the  empirical  evidence  gained  here  indicates 
that  a  graph  with  100  vertices  is  biconnected  99X  of  the 
time  with  6%  edges  whereas  a  graph  with  50  vertices  is 
I  biconnected  99X  of  the  time  when  it  has  \2X  of  the  maximum 

i 

<  vertices  included. 
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RATIO 


5%  10%  15%  20% 

EDGES 

RATIO  OF  LARGEST  BICONNECTED  COMPONENT  TO  ORIGINAL  SIZE 

Figure  17 
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In  investigating  the  effects  of  dividing  on  graph  color¬ 
ing/  note  from  fig.  17  that  with  3G%  edges#  the  largest 
biconnected  component  is  very  close  to  the  original  graph 
for  all  n.  Therefore#  in  order  for  dividing  to  be  benefi¬ 
cial#  random  graphs  must  be  very  sparse.  For  the  empirical 
results  shown  in  fig.  18  graphs  were  generated  with  from  10% 
to  30%  edges.  The  benefit  of  dividing  is  then  shown  by 
plotting  the  ratio  of  the  number  of  expansions  needed  for 
the  divided  graph  to  the  number  of  expansions  for  the  undi¬ 
vided  graph  against  the  percentage  of  edges  present. 
Inasmuch  as  the  relative  benefit  was  the  same  for  all  order¬ 
ing  methods#  only  LCA  was  used  for  this  graph. 

It  appears  as  though  the  benefits  from  dividing  are  in¬ 
significant  since  only  very  sparse  random  graphs  have  arti¬ 
culation  points.  But  many  real-life  problems  which  are 
modeled  by  graphs  may  be  quite  dense  yet  have  articulation 
points.  This  is  true  particularly  in  t ranspor t at i on  and 
communications  problems  where  some  of  the  vertices  may 
represent  junctions  or  hubs  through  which  all  traffic  must 
flow.  Consider  the  city  of  Chicago  in  a  nationwide  tran¬ 
sportation  network  for  example.  In  these  cases  graph  divid¬ 
ing  may  result  in  a  much  faster  colorinq. 

Also  there  are  problems  which  have  graph  models  which 
are  indeed  very  sparse.  Consider#  for  example#  the  exam 
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10  OF  EXPANSIONS  REQUIRED 
DED  GRAPH  /  UNDIVIDED  GRAPH 
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The  vertices  of 


scheduling  problem  in  a  large  university, 
th*  grapn  are  the  classes  and  there  is  a  edqe  between  two 
vertices  if  at  least  one  student  is  enrolled  in  both 
classes.  The  question  then  is:  How  many  exam  periods  are 
needed  to  schedule  all  final  examinations  with  no  conflicts? 
Consider  the  hypothetical  but  realistic  case  of  a  university 
with  200  classes,  and  a  average  of  20  students  per  class. 
Since  many  students  take  the  same  curriculum  together  it  is 
not  inconceivable  to  have  a  graph  with  an  average  degree  of 
each  vertex  of  5.  This  would  result  in  a  graph  of  200  ver¬ 
tices  and  S00  edges  or  2  1/2  X.  This  is  a  very  sparse  graph 
which  i9  likely  to  have  several  articulation  points.  In 
fact,  the  biconnected  components  are  likely  to  be  the  vari¬ 
ous  departments  of  the  university.  Dividing  the  graph  prior 
to  solving  this  problem  would  surely  result  in  a  time  and 
space  savings. 
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IV.  CONCLUSIONS 


The  branch  and  bound  algorithm  described  in  this  paper 
can  be  used  to  find  the  chromatic  number  of  graphs.  Though 
no  empirical  evidence  was  gained  for  graphs  larger  that  52 
vertices#  it  is  felt  that  somewhat  larger  graphs  can  also  be 
colored  in  a  reasonable  amount  of  time  on  the  average  with 
the  LEAST  COLOR  AVAILABLE  method  of  ordering  the  vertices. 
Further#  some  improvement  in  the  lower  bound  function  or 
search  strategy  might  be  possible  which  would  make  the  algo¬ 
rithm  more  efficient.  For  example#  some  lookahead  procedure 
might  conceivably  improve  the  method  of  assigning  costs  to 
each  node  of  the  search  tree.  Also#  the  priority  of  a  node 
in  the  gueue  was  primarily  the  number  of  colors  in  the  par¬ 
tial  coloring  and  only  secondarily  a  function  of  the  depth 
in  the  tree.  A  possible  improvement  is  to  make  the  priority 
some  algebraic  combination  of  the  two  factors  thus  giving 
more  emohasis  on  the  deoth  in  the  tree.  This  might  result 
in  some  faster  colorings.  Since  the  emphasis  here  was  on 
finding  the  best  branching  function#  little  effort  was  ex¬ 
pended  in  looking  for  improvement  in  the  other  2  components 
of  the  branch  and  bound  strategy. 

For  the  five  branching  functions  investigated#  LCA  seems 
far  superior  for  all  graphs  in  both  the  time  reouired  and 
the  amount  of  computer  storage  space  needed.  Though  all 
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five  functions  have  worst  case  time  and  space  complexity 


which  appears  to  be  exponential,  the  LCA  method  has  an  aver¬ 
age  case  time  complexity  which  appears  to  be  nearly  linear 
in  the  range  of  graphs  from  5  to  3 2  vertices. 

Though  the  research  and  experiments  done  here  focused  on 

the  exact  algorithm  for  graph  coloring,  a  few  very  minor 

changes  would  result  in  an  approximation  algorithm.  The 

function  EXPAND  could  be  changed  such  that  it  only  generates 

the  first  possible  child  node  rather  than  all  children  of  a 

node.  No  priority  queue  would  be  needed  and  the  number  of 

nodes  required  in  the  search  tree  would  just  be  n,  the 

number  of  vertices  in  the  graph.  So  the  space-complexity  of 

the  algorithm  would  be  linear  in  n.  The  time  requirea  would 

4 

be  0 ( n  )  if  LCA  were  used  for  an  ordering  function,  but  it 
would  be  interesting  to  see  how  good  its  ability  to  approxi¬ 
mate  the  real  chromatic  number  is  compared  to  other  order- 
i  ngs . 

Dividing  a  graph  at  its  articulation  points  prior  to 
coloring  is  an  effective  way  to  reduce  the  time  required  to 
find  the  chromatic  number.  However,  unless  the  graph  is 
known  to  have  articulation  points  or  unless  it  is  very 
sparse*  the  effort  spent  in  finding  the  biconnected  com¬ 
ponents  may  not  be  worth  it.  *hen  a  graph  does  have  articu¬ 
lation  points,  there  is  a  two-fold  benefit  gained  by  divid- 
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i  ng  the  graoh 


The  first  i$f  of  course*  that  the  size  of 


the  resulting  graph  is  smaller  and  thus  easier  to  color. 
The  second  benefit  comes  from  the  fact  that  the  biconnected 
components  tend  to  be  dense  which  makes  coloring  easier. 


APPENDIX  A 


This  appendix  contains  the  source  listing  tor  the  com¬ 
puter  programs  used  in  this  research.  The  program  is  organ* 
i red  into  various  functions  which  share  some  common  qlobal 
variables.  These  variables  are  listed  in  the  module  EXTERN 
DECS.  The  remainder  of  the  modules  contain  one  or  more 
functions  grouped  together  according  to  their  logical  func* 
tion.  The  modules  and  the  function  in  this  appendix  are: 


MODULE 

EXTERN  DECS 
MAIN 

GENGRAPH 

ORDERNODES 

COLORALL 

DIVIDE 

TREEMANAGER 


FUNCTIONS 

none 
■nainO 
gengra  f ( ) 
random ( ) 
ordernode ( ) 
d  i  c icnoae  ( ) 
col  oral  1  ( ) 
expandt  ) 
c onvb  C ) 
divide!) 
DrocO 
get  node  C ) 
addq  (  ) 
get  t  oo (  ) 
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/****************************** 

GLOBAL  MACRO  DEFINITIONS 

***<***************************/ 
•define  NUMNODES  4000  /*  Maximum  number  of  tree  nodes  */ 

•define  NUMQ  2500  / *  Maximum  size  of  the  priority  Q  * / 

•define  NULL  0  /*  The  empty  pointer  */ 

•define  NUMIRIALS  100  /*  The  number  of  trials  for  each  size 

of  graph  */ 

•define  MAXSIZE  50  /*  The  largest  graph  size  (n)  */ 

•define  M1NSIZE  10  /*  The  smallest  graph  */ 


external  DECLARATIONS 


/*  trnode  is  the  structure  of  the  nodes 
which  make  up  the  search  tree  */ 

struct  trnode  ( 

unsigned  depth  :8; 
unsigned  cost  :8; 
unsigned  vertex  :8; 
unsigned  color  :8» 
struct  trnode  xparptr; 

> ; 

/*  t node  is  a  array  of  search  tree  nodes, 
g  is  an  array  of  pointers  to  tree  nodes  */ 

Struct  trnode  t node  (NUMNODES]  t *q CNUMQ]  } 

/*  ndata  is  a  record  of  data  regarding  the 
graph  vertices  to  be  used  in  ordering  the 
vertices  by  the  branching  function, 
order  is  an  array  of  records  */ 

St  rue  t  ndat  a  ( 

unsigned  next  :  5; 
unsigned  degree  :  5; 
unj i gned  flag  :  1 ; 

>  order  1521; 

int  fd»  /*  file  descriptor  for  output  */ 


i  nt  count ; 
int  hi  node ; 
int  i ndex ; 


int  as i ze » 
int  C  h  rom ; 
int  co 1  or (32) » 
i  n  t  o  r dnum ; 
i nt  saven» 


/*  file  descriptor  for  output  */ 

/*  number  of  vertices  */ 

/*  number  of  calls  to  EXPAND  */ 

/*  first  vertex  to  color.  */ 

/*  index  of  next  available  node  in  the 
array  of  tree  nodes  -  also  the  number 
of  nodes  in  the  search  tree  at  any  time 
/*  index  into  the  priority  queue  */ 

/*  The  chromatic  number  of  the  graph  */ 
/*  The  color  class  of  each  vertex  */ 

/*  The  ordering  function  being  used  */ 
/*  The  saved  value  of  n  for  DIVIDE  */ 
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i n t  graf[32J  [32];  /*  The  adjacency  matrix  */ 

float  nunsuOsi  /*  The  number  of  oiconnected  subgraphs  */ 
float  bigsubn;  /*  The  size  of  the  largest  biconn.  comp.  */ 
float  bigsube  /*  Number  of  edges  in  bigsupn  */ 

long  int  graph [32] ;  /*  The  adjacency  matrix  A  in  bit  strings  */ 

int  sonct [32] ;  /*  A  counter  of  child  nodes  for  each  level  */ 

int  parct[32);  /*  Parent  counter  for  each  level  */ 


/****************************** 

MAIN 

*****************************•/ 
/*  This  function  provides  for  the  flow  of 
control  throughout  the  program  and  records 
the  data  which  is  to  be  collected.  It  also 
contains  the  code  for  any  output  generated  */ 

main  ( )  { 

int  full;  /*  The  percentage  of  edges  */ 

int  edges;  /*  The  number  of  edges  */ 

int  it  it  t*  iteration  variables  * / 

int  ndflaglbi;  /*  Too  large  of  a  search  tree  flag  for  each 

ordering  method  */ 

int  trial;  /*  The  trial  number  for  each  graph  size  */ 

float  maxedges;  /*  The  maximum  edges  possible  for  graph  */ 

float  maxbrlbl#  /*  The  maximum  expansions  reguired  */ 

float  maxndtol#  /*  The  max  nodes  required  */ 

f 1  oat  avgbr  [6]  ;  /*  The  average  expansions  required  */ 

float  avgnd  16] ;  /*  The  average  nodes  in  the  search  tree  */ 

char  *  1 abe 1  [6] ; 

/*  The  output  file  and  labels  */ 
fd  =  c ooen ( " /wo r k / r er au t en/ou t 3 " » ' a  '  )  ; 

1 abe 1  111  =  "RANDOM" ; 

1 abe 1  12]  =  "OEGREE "  ; 

1 abe 1  t3]  =  "DEG  3"; 
label  [«]  =  “ADJCENr; 
label  (5)  =  "LCA" ; 

/*  For  each  value  of  n  desired  */ 
f or (n=MINSIZE;n<sMAXSIZE;n+=5)  < 

/ *  compute  maximum  edges  possible  */ 
maxedges  =  (n  *  (n-i)]/2; 

/*  For  each  value  of  the  percentage  of 
edges  des i red  */ 
for ( full =2; ful 1 <  =  20; f ul 1 +  =  2)  < 


t 
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/*  compute  number  of  edges  */ 
edqes  =  full  *  maxedges/ 1 00 ; 

/*  Print  headings  */ 
print f(fd#"0  =  X3d,  EDGES  =  X4d ( XidXX )  0  , 
n / edges / full); 

printf(fd»"X5sX10s%10sO, "order", "branches"#"nodes"); 
printf(fd,"X15s%12s","avg  max","avg  max  "); 
pr i nt  f ( f d, " 0 ) ; 

/*  Initialize  data  collection  variables  */ 
f or ( i s 1 ; i <s5» i ♦♦ )  t 
maxbr  l  i 1  =  0 ; 
maxnd  l  i )  -  0 ; 
avgb  r  C i 1  =  0 ; 
avgnd  l  i )  =  0 ; 
ndflagti)  =  0; 

> 

/*  For  the  number  of  trials  desireu  / 
for(trial=l;trial<=NUMTRIALS;trial++)  < 

/*  Generate  a  random  graph  */ 
gengraf (n, edges); 

/*  For  each  ordering  method  */ 
for(ordnum=5;ordnum<=5;ordnum+t)  { 
if(i  ndf 1 aq tordnum)  )  < 

/*  Divide  the  graph  and  color  the  components  */ 
i f ( i  divide!))  ndf 1 ag  [ordnumi  =  l; 

/*  Record  data  */ 

a vgb r  lo rdnum)  count; 

avgnd  (o rdnum)  +  =  index; 
ifCcount  >  maxbr  tordnum]  ) 
maxbr tordnum)  =  count; 
iftindex  >  maxnd tordnum]  ) 
maxnd tordnum]  -  index; 

> 

> 

) 

/*  Print  out  the  data  collected  */ 
f or ( i = 1 ; i <=5 ; i *+ )  t 

i  f  ( i  ndflagti))  < 

/*  Compute  the  averages  */ 
avgb r  t i J  /  =  NUMTRIALS; 
avgnd  t i ]  /  =  NUMTRIALS; 
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> 

cex i t ( ) ; 

> 


printf(fd/"X6s"/label  ( i  1  )» 

orintf(fd#"%a.0f%5.Qf"#avgbr(i] »  maxb  r  [  i 1  )  ; 
orintf(fd/"Xb.0f%5.0f  H>avgndtil#maxndlil  ); 
pr i nt  f ( f Oi " 0 ) # 

> 


e  I  se 
1 


pri nt f ( fd, "*bs  NOT  ENOUGH  SPACE  AVAILO, 
label  til  ) ; 


/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

GENERATE  GRAPH 

A*****************************/ 

/*  This  function  generates  a  random  graph  with 
n  vertices  and  e  edges  and  records  it  in  a 
global  adjacency  matrix  */ 

geng raf (n,e)  int  n,e;  < 

i nt  i , j ,a»b; 
long  int  convbt); 

/•  initialize  the  adjacency  matrix  */ 
for  (i=0;i<=n;i+t)  { 

graph  1 1 ]  s  0  ; 

for(j=0;j<sn;jt»)  g  r  a  f  t  i  1  tj  1  -  0 ; 
graf  I  i 1  t i 1  =  1 } 

) 

while  ( e  >  0 )  { 

/*  Get  2  random  vertices  */ 
a  =  r andom ( n  )  J 
b  =  r andom ( n ) i 

/*  Find  the  next  unqenerated  edge  */ 
wh i 1 e (gr a f t a)  tbl  ==  1)  { 

a=(aXn)*i; 
if  (a  ==  1)  b  =  (b  X  n)  ♦  1 
> 

/*  Record  the  new  edge  */ 
graftal  lb]  s  l; 
graftblta)  =  1 ; 
e— ? 

> 

/*  Put  adj  matrix  in  array  of  bit  strings  */ 
for  ( i =1 ; i <sn; i ♦♦)  { 


sa 


graf  Ci 1  l  i  1  =  0; 

for  ( j  =  l  J  j  <sr»;  j  ♦!)  i  f  Cgra  f  t  i  1  [  j  3  )  grapnli]  !=  convb(j); 

> 

> 


/*  ***************************** 

RANDOM 

******************************/ 
/*  This  function  is  used  to  get  random  numoers 
between  1  and  the  input  argument  n. 
taken  from  Grogono  [ref.  --1  */ 

int  random(n)  int  n!  { 

static  long  int  seed  18403; 
int  i  ; 

seed  =  (25173  *  seed  ♦  13849)  X  65536; 
i  =  seed  X  n ; 
ret  urn ( i 1 1 ) ; 

> 


/****************************** 

ORDERNODE 

******************************/ 
/*  This  function  orders  the  vertices  of 
the  graph  for  each  of  the  static  ordering 
methods.  It  does  so  by  filling  the  array 
order. next  with  the  vertex  to  follow  each 
vertex  in  the  graph.  */ 

o  rde  mode  ( )  ( 

int  i r  j . k ;  /*  iteration  variables  */ 

int  save;  / *  the  saved  last  vertex  */ 

int  temphi;  /*  the  temporary  next  vertex  */ 

int  max;  /*  any  maximum  needed  to  maintain  */ 

int  colsum(32J;  /*  used  in  Adjacency  ordering  to  keep 
track  of  number  of  previously  colored 
vertices  a  vertex  is  adjacent  to.  * / 
int  maxcol;  /*  The  column  with  the  max  l's  in  col  sun  */ 

int  hideg;  /*  The  highest  degree  of  any  vertex  */ 

float  dvec  (4)  132);  /*  The  vector  of  degrees  in  DVEC  3  */ 

float  maxdvec;  /*  The  max  vector  calculatea  in  DVEC  3  */ 

long  int  mask;  /*  A  bit  manipulation  mask  */ 
long  int  convbf);  /*  A  converting  function  (SUPP  MOD)  */ 

max  =  0; 
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f  or  ( i 


/*  initialize  the  vertex  data  */ 

-  0 ; i <=n ; i ♦+  )  { 

order t i J . degree  =  0; 
order  ( i ] . nex  t  =  0 ; 
order ( i ] . f 1 ag  =  0 ; 
mask  =  graph  l i ) ; 

/*  Compute  the  degree  of  each  vertex  */ 
f or ( j  =  1 ?  j  <sn> j ♦♦ )  { 

if  ((mask  &  01)  1)  order  li ] .degree  1; 

mask  >>=  1; 

> 


/*  find  the  max  degree  and  the  vertex  */ 
i  f (order (i  1  .degree  >  max)  { 

max  -  orde r ( i 1 . deg r ee » 
h i node  -  i # 

) 

) 


/*  Order  the  vertices  depending  on  the 
current  ordering  number  */ 
swi tch  (ordnum)  { 
case  0: 

break ; 

/*  RANDOM  */ 

case  1: 

f or ( i s0 ; i <n ; i +♦ )  order  (i )  .next  =  i  ♦  1* 
hinode  =  1; 
break ; 

/*  by  OEGREE  */ 

case  2: 

order  [0} .next  =  hinode; 

/*  .flag  indicates  already  colored  vertex  */ 
order (h i node! . f 1 ag  =  1) 
save  -  hinode; 

/*  For  each  vertex  find  the  next  one  */ 
f or ( i s 1 ; i <sn ; i +  ♦  )  i 
max  =  0; 


/*  Find  the  highest  degree  of  the  uncolored 
vertices  */ 


f or ( j  =  1 ; j <  =  n ; j  +  +  )  ( 

i f ( (order l j  1  .degree  >  max)&&(!  order  [  j 1  .  f 1 ag ) ) 


< 


max  =  o rde r  l  j J  . deg r ee ; 
t  emph i  =  j  ; 

> 

> 

/*  record  the  next  vertex  */ 
order  leave]  .next  -  temphi; 
order tt emph i 1 . f 1 ag  =  1 ; 
save  =  temph i f 
> 

break ; 

/*  by  DVEC  3  */ 

case  3: 

/*  initialize  Oth  iteration  to  degree  */ 
f o r ( i = 1 ; i <sn ; i t  +  )  dvecCQlli]  =  order  l  i J  .degree ; 

/*  For  3  iterations  */ 
f or ( i = 1 ; i <®3; i ♦♦ )  { 

/*  For  each  vertex  */ 
for( j  =  ! ; i <=n; j  ♦  +  )  { 
dvec  ( i 1  ( j 1  =  0 ; 

/*  For  every  other  vertex  in  the  graph  */ 
for(k=l ;k<=n?k++)  { 

mask  =  convb(k); 

/*  If  they  are  connected  compute  summation 
for  the  ith  iteration  */ 

if((graph(j]  &  mask)  >  0) 

dvecCiJtjJ  +=  dvec  (i'll  [k]; 

> 

> 

) 

/*  Order  the  vertices  */ 

/*  The  ordering  method  is  identical  to  case  2 
except  for  the  ordering  parameter  */ 
save  =  0; 

f  or ( i s i ; i <*n ? i ♦ t )  { 
maxdvec  =  0; 
f  or ( j  =  l ; j  <sn ; j *  *  )  < 

i f ( (dvec f 3)  [  i  J  >  maxdvec)  &&  (1  order  ( j ). f 1 aq ) )  { 
maxdvec  =  dvecl3)  [jj? 
t  emph i  =  j ; 

) 

} 
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order ( save! . ne x t  =  temphi; 
order (temph i 1 . f 1 ag  =  l; 
save  =  t  emph i ; 

> 

hinode  =  order (0) ,ne«t ; 
break; 


/*  by  ADJACENCY  */ 

case  4: 

/*  initialize  */ 
f or ( i si ; i <*n; i ♦♦ )  colsumti) 
order  10] . nex t  =  Hinode; 
order th i node] . f 1 ag  =  l; 
save  -  hinode; 
col  sum  (save)  -  =  33; 


o; 


/*  For  each  vertex  find  the  next  one  */ 

f or ( j *  1 ; j <sn-l ; j ♦♦)  ( 

max  =  0; 
mask  =  0 1 ; 

/*  increment  the  bit  string  col  Sum  with  each 
vertex  added  to  the  set  of  colored  vertices  »/ 
f or ( i  s 1 ;  i  <*n ;  i  ♦♦ )  ( 

i f ( (graph  [save]  &  mask)  >  0) 
col  sum [i ] 


i  f  (col  sum  (i  J  >  max)max  =  colsumli]  ,* 
mask  <<=  l; 

> 


/*  Order  the  vertices  -  again  the  same 
ordering  method  is  used  as  in  case  2  except 
the  ordering  parameter  is  by  value  of  col  sum  */ 
hi  deg  =  0; 
f o  r ( i = 1 ; i <sn ; i ♦♦ )  ( 

i f ( (co 1  sum [ i  1  ==max)&&(orderli]  ,degree>hideg) 

&i(i  o rde r ( i ] . f 1 ag ) )  i 
h i deg  =  order  (i ] .degree; 
t  emph i  s  i  ; 

> 

} 

order  (save)  . ne x t  =  temphi; 

order [t emph i ]. f 1 ag  =  i; 

col  sum  (temoh i ]  -=  33; 

save  -  t emph i ; 

) 

break; 
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> 


/*  by  LEAST  COLORS  AVAILABLE  */ 

case  5: 

/*  no  work  done  in  this  function. 
See  PICKNOOE  */ 
break ; 


/*****************************» 

PICKNODE 

****************************** 
This  function  is  called  from  the  routime 
EXPAND  to  get  the  next  vertex  to  color. 

If  one  of  the  static  orderings  is  being 
used  it  just  looks  uo  the  next  vertex  in 
the  table  order. next.  If  the  dynamic 
ordering  LCA  is  being  used,  it  must  trace 
the  search  tree  to  the  root  to  determine 
the  next  vertex  to  color.  */ 

p i c knode ( nodept r )  struct  trnode  »nodeotr;  i 

/*  Declarations  */ 
long  int  cc 1 32] , vdone , 1 ong 1 ; 
int  max,isum,k,i,sum(321 »nodenum,hiaeg; 
struct  trnode  »tempotr; 

/ *  If  static  ordering,  just  table  lookup  */ 
ifCordnum  <  5)  return (order (nodept r»>vertexj .next); 

/ *  initialize 

vdone  is  a  bit  string  to  record  which 
vertices  have  already  been  colored, 
cc  is  the  color  class  matrix 
sum  counts  the  colors  not  avail  for  eacn 
vertex.  */ 

vdone  =  0, 
f or ( i - 1 ? i <  =  n  ;  i ft )  { 
cc ( i 1  -  Or 
sum  l  i ]  =0  # 

> 

tempptr  =  nodeotr; 

/*  Trace  from  node  N  to  root.  */ 
while  (tempptr  NULL)  { 

/*  fill  color  class  matrix  */ 
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cc ftemppt r->col or)  !=  graph  [ t emppt r-> ve r t e *  1 ; 

/*  record  colored  vertices  */ 
vdone  !=  convb ( temppt r->ver te* ) ; 
tempptr  =  t emppt r->parpt r ; 

} 

1 ong 1  =  015 

max  s  0# 


/•  For  each  of  the  n  vertices  */ 

f or (k  =  l ; k<=n?k*t)  < 

/*  If  the  vertex  k  is  not  colorea  */ 
if((vdone  &  convb(k)  )  0)  { 

i  9 urn  s  0 ; 

/*  count  the  color  classes  to  which  vertex 
k  cannot  be  added  */ 
for ( i *1 ; i <snodept r->cost ; i ♦  ♦  )  < 

isum  t s  longl  &  (cc(i)>>(k~l)); 

> 

Sum(k)  =  isum; 

/*  keep  track  of  fewest  number  of  color 
c 1  asses  avai 1 abl e  */ 
ifCiSum  >  max)  max  s  isum/ 

> 

> 

hideg  =  0; 

/*  Of  all  the  vertices  with  the  fewest  color 
classes  available*  find  the  one  with  the 
highest  degree  */ 
f or ( i =1 ; i <=n; i tt )  < 

i f ( ( sum  [  i ]  ==  max)  &&  (J(convbCi)  &  vdone)))  < 

i f (order  l  i 1  .degree  >  hideg)  < 
hideg  s  o rder  ( i ]. degree f 
nodenum  -  i } 

) 

> 

> 

/ *  Return  the  vertex  number  of  the  next 
vertex  to  color.  */ 


return(nodenum) * 
) 


/*  ***************************** 

COLOR  ALL 

****************************** 
This  function  creates  the  root  noae  of  the 
search  tree,  outs  it  in  the  oriority  queue 
then  calls  the  function  EXPAND  to  exoana 
the  top  node  in  the  queue  until  the  top  of 
the  queue  is  a  complete  coloring.  */ 

co 1  ora  I  1 ( )  { 

struct  trnode  *nodept r , *get t oo ( ) ? 
i  nt  i  ; 


/*  Initialize  for  new  coloring  */ 

k  =  o; 

index  s  0; 
qs i ze  -  Or 
count  =  Or 
tnode  101 .cost  =  31# 

for  ( i  =0;  i  <N(JMQ;  i  *♦)  qCil  =  &tnode[0]; 


/*  Get  an  empty  node*  fill  it  with  the 
data  for  the  root  of  the  tree  and 
out  it  in  the  queue.  */ 
if  Cl  Cnodeptr  =  getnodeC)))  return(O); 
nodept r->parot r  *  NULL; 
nodeot r->deot h  s  l; 
nodept r*>cost  a  1; 
nodept r->vert ex  =  hinode* 
nodeot r->co I  or  =  It 
addqCnodept  r ) ; 

/*  ' forever  *  */ 

for(;;i  { 


/*  Get  the  too  node  in  the  queue  */ 
nodeotr  =  qettop()» 


/*  If  this  is  complete  coloring,  quit  */ 
i f (nodeot r->Oeot h  ~~  n)  break? 
else  if  (1  exDand (nodeot r  )  )  return(O); 

> 


/*  Record  the  chromatic  numper  */ 
chro#  -  nodept r->cost ? 

/*  Record  the  color  class  assignment  of 
each  vertex  in  this  Dath.  */ 
while  (nodept r->oarot r  1=  NULL)  ( 

col  or [nodeot r->vertexl  -  nodeot r->co 1  or ; 
nodeotr  s  nodeot r->oarot r ; 
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> 

co 1  or (nodeot r«> ver t ex ]  =  nodept r->co 1  or ; 

return(l); 

> 


/****************************** 

EXPAND 

****************************** 

This  function  expands  the  node  pointed  to  by 
the  variable  nodeptr.  It  creates  a  child 
node  for  each  possible  color  available  as 
well  as  a  node  for  one  new  color.  */ 

expand  (nodeptr)  struct  trnode  *nodeptr;  { 

struct  trnode  * t emppt r , ‘newnode * *ge t node ( ) ; 

i nt  vertexnum, i } 

long  int  mask , co 1 ava i 1 , con vb ( ) ; 

/*  increment  counters  */ 

count  *-  1; 

caret  (nodept  r->depthi  *■*■! 

/*  Get  the  number  of  the  next  vertex  to  color  */ 
vertexnum  =  o i c knode ( nodeP t r ) ; 

/ *  create  new  nodes  for  all  colorings  using 
only  existing  colors  */ 

co 1 ava i 1  -  01; 

for  ( i - 1 ; i <nodept r->cost ; i ♦  ♦  )  ( 

col avai 1  <<  =  1  ; 
co 1 ava i 1  I  -  0  1  } 

/*  colavai!  is  a  bit  string  with  a  one  in  the 
position  of  every  available  color.  */ 

> 

tempotr  =  nodeotr; 

/*  For  each  node  N  i i  the  oath  to  the  root  */ 

for  ( ; ; )  { 

/*  If  N. vertex  and  vertexnum  are  connected  */ 
i f ( (graph  (vertexnum)  4  c on vb ( t emppt r->ve r t e x  )  )  1=  0)  ( 

mask  =  convb( temppt r->co 1  or  ) ; 

/*  Remove  N. color  from  col  aval  1  */ 
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H 


if  ((colavail  i  mask)  1  = 
/ 

if  ( t emoot r->parpt r  =  =  NULL) 
tempptr  =  t empp t r->pa ro t r ; 

> 

/ *  If  any  colors 
if  (colavail  is  0)  { 
mask  =  01; 


0 )  colavail 
break; 


are  left  * / 


t  s  mask ; 


/*  For  each  color  available  */ 
for  ( i =1 ; i <=nodept r->cost ; i +t)  ( 
if  ((colavail  4  mask)  is  Q)  ( 


/*  Generate  child  node  */ 
if  l (newnode  =  getnoaet ) )  =  =  0)  return(O); 
newnode->deDt h  =  nodept r->deot h  t  l; 
newnoae->cost  =  nodect r->cost ; 
newnode*>verte»  =  vertexnum; 
newnode*>co 1  or  =  i; 
newnode*>oarctr  s  nodeotr; 

/*  increment  child  counter  */ 
sonct  lnodeotr->depth] 


/*  Put  the  node  in  the  queue  */ 
addq ( newnode ) • 

) 

mask  s<<  1 ; 

) 


} 

/*  Generate  one  more  node  with  a  new  color 
if  ((newnode  s  getnodeO)  =  =  0)  return(O); 
ne«noce*>dfDt h  =  nodept r->deo t h  +  l; 
newnode*>cost  =  nodept r->cost  t  t; 
newnode->verte*  -  vertexnum# 
newnode->co I  or  -  nodeot r->cost  *  l; 
newnode->parDt r  =  nodeotr) 
sonct  (nodeot r->depth) 
addq (newnode ) ; 
ret  urn  (  1  ) ; 


> 


*/ 


/*  A**************************** 

CONVERT  TO  BIT 

a****************************  */ 

f*  This  function  converts  an  integer  argument 
between  0  and  16  to  a  bit  string  which  has  a 
1  in  the  position  of  the  argument  and  0's 
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everywhere  else  */ 
long  int  convb(ary)  int  arg;  { 

»  n  t  i  i 

1 ong  int  b i t  oos  J 
bi tPOS  =  0 1 ; 

f or ( i = 1 ; i <arg; i )  bitpos  =  <<  1; 
returnCbi  toos)  ; 

> 


/  * 


*********************** 

DIVIDE 

**********************/ 
/*  This  function  divides  the  graph  into  its 
biconnected  components.  It  then  creates  a 
new  adjacency  matrix  f«r  each  component  and 
calls  the  function  COLQRALL  to  color  the 
component.  It  records  the  data  for  each 
component  and  sums  it  up  for  all  comps  of  a 
graph.  When  it  returns*  the  sum  data  is  in 
the  variables  set  up  by  MAIN.  */ 


di  v  i 

de()  ( 

i  nt 

number  C 32} ; 

/  * 

i  nt 

back  132]  ; 

/* 

i  nt 

stackv  (32] ; 

/* 

i  nt 

num  * 

/* 

i  nt 

W(VitOP? 

/* 

i  nt 

st  vpt  r ; 

/  * 

i  n  r 

st  eot  r ; 

/* 

i  nt 

nume* 

/* 

i  nt 

num  v  } 

/* 

i  nt 

vnum [32]  } 

/* 

t 

i  nt 

totcount  ; 

/  * 

i  nt 

t  ot i ndex ; 

/* 

i  nt 

t  ab 1 e  (321  ; 

/  * 

i  nt 

t  ab 1 ept  r ; 

/* 

i  nt 

vert l *  vert  2  * 

/  * 

1  ong 

i  nt  convb  (  )  ’• 

i  nt 

i  *  j  J 

of  vertices*/ 
hich  a  oack 


to  record  the  visiting  order 
record  the  lowest  vertex  to 
edge  exist  */ 

A  stack  of  the  vertices  */ 

To  maintain  the  next  number  */ 
vertices  -  too  is  too  of  stack  */ 

Pointer  to  top  of  vertex  stacx  */ 
oointer  to  too  of  edge  stack  */ 
number  of  edges  in  bi-component  */ 
number  of  vertices  in  b i -c omoonen t 
array  to  record  the  vertex  mumpers 
he  edges  of  the  edge  stack  are  popped  */ 
count  the  total  calls  to  ExPAND  for  all 
components  */ 

count  the  total  nodes  used  */ 
to  convert  the  bicomoonent  to  a  graph  */ 
Pointer  into  table  */ 

the  vertices  of  the  edges  in  new  graph  */ 


/*  An  edge  in  a  biconnected  component  defined 
by  the  two  vertices  at  its  end  */ 


*/ 

a  s 
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struct  edge  < 

uns i gned  v 1 :  8 ; 
unsigned  v2:  8; 

>  s t  ac ke  1 3501  ; 

/*  Stacke  is  an  stack  of  edges  */ 

/*  Save  off  the  oringginal  value  of  n  -  the 
biconnected  components  will  be  different  »/ 

seven  ~  n» 

/*  Initialize  all  variables  */ 

number (11  =  1 ; 

num  =  2 ! 
s  t  ac  k  v ( 01  =  1 ; 

too  ~  1  * 

f o r ( v  =  2 } v<  =  n ; v ♦  +  )  number  (vl  =  0; 

steot  r  =  -1 ; 
stvptr  s  0; 
totcount  =  0; 
tot i ndex  =  0 ; 

/*  Begin  a  deoth  first  search  of  the  graph 
exploring  forward  edges  */ 

forward : 

/*  While  there  is  an  unexplored  forward  edqe 
from  the  top  vertex  on  the  stack  to  some 
vertex  w  */ 

while((w  =  proc(topl)  0)  < 

/*  Stack  the  edge  (toprw)  on  the  edge  stack  */ 
stackeft+steptr] .vl  =  top; 
stacke  (stept rl  . v2  =  w# 

/ *  If  w  has  already  been  explored  give  the  min 
of  number(w)  and  back(top)  to  bac<(toc). 
i f ( number  (wl  >  0  )  { 

if(number(wl  <  back(tooj)  b  a  c  k  ( t  o  o  1  =  n  u  m  d  e  r  f  w  )  ; 

> 

/*  Else  stack  w  and  explore  it  */ 

else  i 

number twl  s  num+t; 
bac  k (wl  ~  numoertwi; 
stackv  tffstvotr]  -  w; 
top  =  w; 

> 

1 

/»  If  there  is  no  forward  edge  then  bacxuo 
to  the  last  vertex  which  still  has  a  fcrwaro 


95 


/  AO-A097  330  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA  F/G  9/2 

AN  INVESTIGATION  OF  SEVERAL  BRANCHING  FUNCTIONS  IN  A  BRANCH  AND— ETC (U) 
DEC  80  R  E  RAUTENBER6 

UNCLASSIFIED  NL 


edge  which  is  unexplored  */ 


backup: 

i f ( st vpt r  >  0)  { 

/*  try  the  next  to  top  vertex  */ 
v  -  stackv  1st vpt r  -  1J  » 

/*  if  it  has  a  number  less  that  the  back  of 
the  too  vertex  then  we  have  an  articulation 
point  and  the  bi  connected  component  is  on 
top  of  the  edge  stack  */ 
if(backCtop)  >*  number(vl)  < 

/*  increment  and  initialize  counters  */ 
numsubs  l; 
nume  =  0# 
numv  =  Or 

for( i=l ? i<ssaven; i ♦♦)  graphti]  s  0; 
tableptr  -  0# 
tabletO)  =  0; 

for C i =0? i <ssavenj i ++)  vnum[i]  -  0; 

/*  poo  edges  until  end  of  component  */ 
wh i 1 e( (number (stac ke  (stept rj  .  v  1  j  >*  number(vl)  && 
(number Istacke Isteot r]  .v21  >=  number  tv) ) 

&&  (steotr  >s  0))  ( 

t*  increment  data  counters  */ 
numet* » 

i f ( ttvnum  1st ac ke  (st eot r ]  . v 1 1  ss  l)  numv  +  +  » 
i f ( t  +  vnum  (stacke  (steot r)  ,v2]  ==  1)  numv+t  » 
vertl  s  0J 
vert2  s  0; 

/*  check  to  see  if  either  vertex  is  already 
in  the  table  -  if  so#  get  the  index  which 
is  going  to  be  the  new  vertex  number  */ 
f  or  ( i  =0 » i  <3t abl ept r ; i ++ )  ( 

if(tabie(i)  --  stacke (stept rl . v  1 )  vertl  3  j; 
if(table(il  =s  stacke (steot r] . v2)  vert2  s  i; 
} 

/*  if  not  in  table#  put  it  there  and  get  new 
vertex  number  */ 
i f ( vert  1  ==  0)  ( 
tableptr+t? 

t ab l e (t abl eot rl  =  stacke  (stept r] .vl # 
vertl  3  tableotr# 

> 

i f ( vert2  ==  0)  { 
tableptrtr# 

tabl e  (t abl ept rl  3  stacke (stept rl ,v2# 
vert2  s  tableptr; 
l 
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/*  record  the  edge  between  the  new  vertex 
numbers  */ 

graphlvertll  }a  convb(vert2) ; 
graph(vert21  5*  convb ( vert  1 ) t 

y*  decrement  edge  stack  pointer  and  do  again  */ 
stept  r— j 
> 

/*  The  biconnected  component  is  now  recorded 
in  the  adjacency  matrix  and  ready  to  color  */ 

/*  Is  this  the  biggest  component  of  this  graph  */ 
if(numv  >  bigsubn)  bigsubn  a  numv; 
if(nume  >  bigsube)  bigsube  •  numef 

/*  new  n  is  site  of  the  table  */ 
n  s  tableptr? 

/*  Call  the  ordering  function  */ 
ordernodeO  ? 

/*  Color  the  comoonent  */ 
if(J  colorallO)  retum(O); 

/*  Increment  the  cummulative  counters  */ 
totcount  ♦=  count; 
totindex  ta  index# 

> 

/*  If  no  art.  pt.  pop  the  vertex  stack  and 
explore  any  forward  edges  of  the  vertex 
which  is  now  on  top  */ 

else  if(back(topl  <  back  (vl  )  backCv)  =  backltop); 
top  a  stackv  (  —  at  vot  rl  ; 
goto  forward; 

> 


/*  When  all  the  connected  vertices  have 
been  explored*  check  for  another  unconnected 
component  and  explore  it  */ 
f or C i =2; i <=saven; i ♦> )  { 
i f (number ( i J  aa  0 )  ( 
numberlil  a  l; 
num  a  2 } 
stackv  CO)  a  jf 
too  a  i; 
stvptr  a  0; 
steptr  =  -t; 
goto  forward; 

> 

> 

/*  When  everything  is  done  and  all  components 
are  found  and  colored*  set  the  counters  to  the 
cummulative  totals*  reset  n»  and  reset  the 
adjacency  matrix  */ 


count  *  totcount; 

Index  =  totindex; 
n  =  seven# 

for  ( i =1 ; i <=n; i  ♦♦)  ( 
grafli)  C 1 ]  =  0; 
for  ( j S1 » j <=n# j ♦♦)  { 

if  (grafli]  Ij)  ss  2)  grafli]  t j  1  *  l; 
if  CgrafCilCjl  =s  1)  grapMi)  {=  convb(j); 
> 

> 

return  Cl); 

> 


/*»***» 

******* 


PROC 

***********************/ 
/*  This  function  determines  if  we  can 
proceed  in  the  depth  first  search  exploring 
forward  edges.  It  marks  the  adjacency 
matrix  with  a  2  as  each  edge  is  explored 
so  that  no  edge  is  explored  twice  */ 


proc(i)  int  ♦;  ( 


int  j; 

/*  Is  there  a  vertex  adjacent  to  vertex  i 
on  an  edge  which  hasn't  been  explored.  If  so 
mark  it  explored  and  return  the  vertex  number  */ 
for(jsi;j<ssaven;j*+) 

i  f  (graf  UJ  Cjl  *=  1)  C 
graf  li ]  ( j  J  =  2; 
graf ( j ) (i)  *  2; 
returnC  j ) # 

> 

return(O) ; 

} 


/ 


TREE  MANAGER 


/ 


/*  This  module  contains  the  functions 
necessary  to  manage  the  search  tree  nodes 
as  well  as  the  priority  queue.  */ 


/*  *********  * 


********** 

GETNODE 


********* 
*******  */ 


/*  This  functions  maintains  the  array  of 
nodes  and  returns  a  pointer  to  the  next 
available  node  when  called.  Index  is  a 
global  variable  */ 
struct  trnode  *getnode()  < 

/*  index  gets  initialized  in  COLORALL  */ 

index  1* 

if  (index  <  NUMNODES)  return(4tnode  (index! ) ; 
else  return(O); 

> 


/*  ***************************** 

ADD  TO  QUEUE 

******************************/ 
/*  This  functions  adds  a  new  node  to  the 
priority  queue.  It  puts  it  in  the  last 
position  and  then  'sifts  it  up'  until  it  is 
in  its  proper  location  */ 
addq  (otr)  struct  trnode  *ptr?  ( 

struct  trnode  *temp; 
i  n  t  i  0 }  ; 

/*  Put  the  new  node  at  the  bottom  */ 

qsize  ts  l; 

if  (qsize  **  NUMQ)  pr i nt f ( f d, "Q  0VERFL0W0); 
qlqsizel  =  ptr; 

/*  j  is  the  node  being  sifted  up  */ 

i  =  qsize; 

/*  i  is  the  parent  of  )  */ 

t  s  j /a; 

/*  While  j  is  not  the  root  */ 

whi 1 e( J  >  1)  ( 

/*  compare  priority  of  j  to  i  */ 
if(qlil->cost  <  q(jJ->cost)  return; 
else  if  (qlil->cost  =  =  q(jJ->cost)  ( 

if  (q(il->deoth  >=  qljl->depth)  return; 

> 

/*  Exchange  nodes  and  continue  */ 

t  emp  =  q ( j 1 ; 
q  C  j  1  a  q  ( i  1  ; 
q  1  i 1  a  temp; 
i  *  i; 
i  *  j/2 f 
) 

> 
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GET  TOP 

******************#***********/ 

/*  This  functions  returns  a  pointer  to  the 
tree  node  which  is  on  top  of  the  priority 
queue.  The  queue  is  a  heap.  When  the  root 
of  the  fuff  binary  tree  is  removed#  the  node 
in  the  last  position  is  moved  to  the  root  and 
then  'sifted  down*  */ 
struct  trnode  *gettop()  < 


struct  trnode  *savept r, *temo) 
i  nt  i #  i # 

/*  save  off  the  root  node  */ 

saveptr  -  q(li; 

/*  move  the  last  node  to  the  root  */ 

q(l]  s  qfqsize)  ? 

/*  decrement  the  qsize  */ 

qsize  =-  l; 

/*  sift  the  root  node  down  */ 

/*  i  is  the  parent  node  being  looked  at  */ 

i  =  1) 

/*  j  is  the  left  child  of  i  */ 

j  =  2  *  i? 

/*  While  were  not  out  of  the  heap  */ 
whi  1  e( j  <*  qsi  ze)  ( 
if  <>  <  qsize)  < 

/*  Get  the  child  of  i  with  least  cost  */ 
if  (qIjl->cost  >  q  [  j +  1) »>cost)  j  =  j*l) 

/*  if  both  the  same  cost#  get  greater  depth  */ 
else  if  (q[jl->cost  *=  q  ( j ♦ 1 1 ->cost ) 

i f (q (j 1 ->depth  «  q ( j  + 1] •>dept h )  J  =  j*l) 

> 

/*  If  both  children  have  lower  priority#  quit  */ 
if  (qCi]->eost  <  q[jl*>cost)  break) 
else  if  (q(il->cost  as  qIjl»>cost) 

if  (q li 1 *>deoth  >s  q(jl*>depth)  break) 

/*  exchange  node  i  with  JCits  highest 
priority  child)  */ 


temp  s  q{ j] ; 
q  l  J 1  3  q  ( U  ; 
q(i)  s  temp) 


f*  set  i  *  to  J  and  do  it  again  */ 


i  s  2  *  i; 

> 

/*  when  done#  return  the  save  off  root  */ 

return Csaveot  r ) ; 

> 


! 

I 


101 


APPENDIX  B 


This  appendix  contains  output  listings  front  the  various 
programs  run  to  generate  the  plots  described  in  Section  IX. 
The  first  listing  is  the  data  for  graph  coloring  without  di¬ 
viding  for  graphs  with  from  10  to  30  vertices.  The  edge 
velues  are  30%»  50X?  and  70%. 

The  second  listing  is  the  output  from  the  DIVIDE  routine 
which  divides  the  graphs  and  records  data  regarding  the  size 
of  the  biconnected  components  and  the  number  of  components. 
The  graphs  ranged  in  size  from  6  to  100  vertices. 

The  next  listing  contains  data  from  tne  routine  to 
divide  the  graphs  and  color  the  components.  Graphs  with  25 
vertices  and  2%  to  1«%  edges  were  used  in  this  experiment. 

The  last  listing  shows  the  effects  of  the  amount  of 
edges  when  n  is  held  constant.  For  this  experiment  n  s  20 
while  the  edges  vary  from  20%  to  BOX  by  5X  intervals. 
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****************************** 
COMPARISON  OF  5  ORDERING  METHODS  FOR  SEVERAL  GRAPH  SIZES 
****************************** 


N  =  10,  EOGES  s  13(30%) 


ORDER 

BRANCHES 

NODES 

avg 

max 

avg 

max 

RANDOM 

19 

65 

38 

115 

DEGREE 

10 

14 

21 

29 

DVEC  3 

9 

11 

21 

25 

ADJCEN 

9 

9 

21 

24 

LC  A 

9 

9 

21 

24 

BRANCH 

FACTOR  BY 

LEVEL 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RANDOM 

1.55 

1  .92 

1.96 

2.14 

1.75 

2.11 

1 .85 

2.43 

2.01 

0.00 

DEGREE 

1.23 

1.59 

1.82 

1  .94 

2.29 

2.41 

2.57 

2.78 

2.88 

0.00 

DVEC  3 

1.04 

1.45 

1.96 

2.22 

2.45 

2.47 

2.70 

2.83 

2.87 

0.00 

ADJCEN 

1.00 

1.29 

2.16 

2.44 

2.41 

2.45 

2.68 

2.83 

2.89 

0.00 

LC  A 

1.00 

1.29 

2.15 

2.29 

2.47 

2.61 

2.74 

2.87 

2.92 

0.00 

N  =  10, 

EDGES  = 

22(50%) 

ORDER 

BRANCHES 

NODES 

avg 

max 

avg 

max 

RANDOM 

29 

129 

53 

239 

DEGREE 

10 

18 

21 

35 

DVEC  3 

9 

11 

20 

25 

ADJCEN 

9 

1 1 

20 

25 

LCA 

9 

9 

20 

25 

BRANCH 

FACTOR  BY 

LEVEL 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RANDOM 

1.34 

1.81 

1.72 

1.99 

1.63 

2.24 

1.57 

2.13 

1.75 

0.00 

DEGREE 

1.08 

1.36 

1.59 

1.76 

1.93 

2.25 

2.57 

2.66 

3.13 

0.00 

DVEC  3 

1.01 

1.18 

1.42 

1.97 

2.27 

2.42 

2.56 

2.73 

3.13 

0.00 

ADJCEN 

1.00 

1.00 

1.35 

2.13 

2.44 

2.50 

2.65 

2.80 

3.16 

0.00 

LCA 

1.00 

1.00 

1.35 

2.07 

2.35 

2.45 

2.72 

2.87 

3.24 

0.00 

N  s  10, 

EDGES  = 

31 (70%) 

ORDER 

BRANCHES 

NODES 

avg 

max 

avg 

max 

RANDOM 

33 

160 

56 

250 

DEGREE 

10 

17 

19 

27 

DVEC  3 

9 

16 

18 

26 

ADJCEN 

9 

10 

18 

22 
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BRANCH  FACTOR  BY  LEVEL 


l 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RANDOM 

1.08 

1.50 

1.48 

1.98 

1.51 

2.07 

1.43 

2.33 

1.48 

0.00 

DEGREE 

1.01 

1.12 

1.24 

1.36 

1.58 

1.98 

2.31 

2.63 

3.10 

0.00 

DVEC  3 

1.02 

1.05 

1.19 

1.21 

1.61 

2.15 

2.51 

2.73 

3.14 

0.00 

ADJCEN 

1.00 

1.00 

1.00 

1.13 

1.78 

2.41 

2.53 

2.72 

3.14 

0.00 

LCA 

1.00 

1.00 

1.00 

1.13 

1.78 

2.23 

2.52 

2.76 

3.25 

0.00 

N  =  15, 

EDGES  = 

31(30%) 

ORDER 

BRANCHES 

NODES 

avg  max 

avg  max 

RANDOM 

NOT  ENOUGH  SPACE  AVAIL 

DEGREE 

22 

88 

48  173 

DVEC  3 

1  7 

42 

38  82 

ADJCEN 

15 

34 

36  66 

LCA 

14 

22 

35  44 

BRANCH  FACTOR  BY 

LEVEL 

1 

2 

3  4 

5 

6 

7 

8 

9 

10 

RANDOM 

NOT 

CALCULATED 

DEGREE 

1.21  1 

L  .59 

1.67  1.86 

2.03 

2.03 

2.09 

2.22 

2.36 

2.49 

DVEC  3 

1.06  1 

..24 

1.77  2.07 

2.11 

2.28 

2.31 

2.35 

2.53 

2.73 

ADJCEN 

1.00  1.02 

1.64  2.13 

2.45 

2.48 

2.47 

2.46 

2.44 

2.65 

LCA 

1.00  1 

1.02 

1.62  2.07 

2.36 

2.41 

2.44 

2.57 

2.59 

2.71 

11  12 

13  14 

15 

RANDOM 

NOT 

CALCULATED 

DEGREE 

2.66  2.83 

3.05  3.44 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

DVEC  3 

2.77  2.88 

3.20  3.43 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

ADJCEN 

2.78  2.98 

3.15  3.44 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

LCA 

2.82  ; 

5.07 

3.32  3,60 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

N  s  15,  EDGES  s  52(50%) 


ORDER  BRANCHES  NODES 


avg 

max 

avg 

max 

RANDOM 

NOT 

ENOUGH 

SPACE  A 

DEGREE 

24 

72 

47 

129 

OVEC  3 

19 

47 

38 

89 

ADJCEN 

16 

61 

35 

115 

LCA 

15 

33 

32 

66 

BRANCH 

FACTOR  BY 

LEVEL 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RANOOM 

NOT 

CALCULATED 

DEGREE 

l  •  l  a 

1.30 

1.47 

1.57 

1.72 

1.81 

1.85 

2.07 

2.30 

2.25 

DVEC  3 

1 .04 

1.17 

1.34 

1.59 

1.98 

l  ,9b 

2.04 

2.21 

2.22 

2.45 

ADJCEN 

1.00 

1.00 

1.04 

1.49 

2.02 

2.43 

2.39 

2.41 

2.25 

2.41 

LC  A 

1.00 

1.00 

1.04 

1.47 

1.95 

2.21 

2.21 

2.29 

2.37 

2.50 

11 

12 

13 

14 

15 

RANDOM 

NOT 

CALCULATED 

DEGREE 

2.32 

2.57 

2.71 

3.14 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

DVEC  3 

2. 49 

2.55 

2.70 

3.16 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

ADJCEN 

2.53 

2.63 

2.87 

3.22 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

LC  A 

2.5a 

2.70 

2.95 

3.38 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

N  =  15, 

EDGES  = 

73(70*) 

ORDER 

BRANCHES 

NODES 

avg  max 

avg  max 

RANDOM 

NOT  ENOUGH  SPACE  AVAIL 

DEGREE 

24 

75 

44  128 

DVEC  3 

21 

56 

38  96 

ADJCEN 

15 

32 

31  66 

LC  A 

14 

28 

29  55 

BRANCH  FACTOR  BY 

LEVEL 

1 

2 

3  4 

5 

6 

7 

8 

9 

10 

RANDOM 

NOT 

CALCULATED 

DEGREE 

1.02  1 

1.07 

1.15  1.31 

1.42 

1  .60 

1.71 

1.78 

1.98 

2.20 

DVEC  3 

1,02  1 

1.02 

1.06  1.24 

1.36 

l  .60 

1.77 

2.01 

2.05 

2.33 

ADJCEN 

1.00  1 

l  .00 

1.00  1.00 

1,08 

1.49 

2.03 

2.31 

2.40 

2.49 

LC  A 

1.00  1 

1.00 

1.00  1.00 

1.08 

1.47 

1.93 

2.14 

2.14 

2.35 

11  12 

13  14 

15 

RANDOM 

NOT 

CALCULATED 

DEGREE 

2.21  < 

>.52 

2.70  3.18 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

DVEC  3 

2.27  2.46 

2.68  3.15 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

ADJCEN 

2.45  2.54 

2.85  3.14 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

LCA 

2.49  2.63 

2.9!  3.44 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

N  s  20, 

EOGES  = 

57(30*) 

ORDER  BRANCHES 
avg  max 


NODES 
avg  max 


mrnrnmm OPPW 


RANOOM 

NOT 

ENOUGH  SPACE  AVAIL 

DEGREE 

63 

312 

125 

601 

DVEC  3 

33 

118 

70 

229 

AOJCEN 

25 

208 

57 

418 

LCA 

22 

121 

53 

247 

BRANCH 

FACTOR  BY 

LEVEL 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RANDOM 

NOT 

CALCULATED 

DEGREE 

1.53 

1  .59 

1.72 

1  .82 

1.86 

1.60 

1.76 

2.08 

2.08 

2.  1  1 

DVEC  3 

1.25 

1.34 

1.63 

1 .87 

2.01 

2.05 

2.06 

2.28 

2.41 

2.32 

ADJCEN 

1.00 

1.00 

1.49 

1.84 

2.23 

2.22 

2.42 

2.40 

2.70 

2.7o 

LCA 

1.00 

1.00 

1.49 

1.77 

2.07 

2.21 

2.36 

2.45 

2.51 

2.59 

11 

12 

13 

14 

IS 

16 

17 

18 

19 

20 

RANDOM 

NOT 

CALCULATED 

DEGREE 

2.17 

2.30 

2.41 

2.49 

2.51 

2.56 

2.76 

2.80 

3.39 

0.00 

DVEC  3 

2.36 

2.37 

2.38 

2.51 

2.57 

2.68 

2.79 

2.91 

3.26 

0.00 

ADJCEN 

2.70 

2.54 

2.38 

2.50 

2.55 

2.59 

2.81 

2.86 

3.31 

0.00 

LCA 

2.52 

2.54 

2.50 

2.56 

2.65 

2.78 

2.98 

3.18 

3.50 

0.00 

N  s  20, 

EDGES  s 

95(50X) 

ORDER  BRANCHES 

NODES 

avg 

max 

avg 

max 

RANDOM 

NOT 

ENOUGH  SPACE  AVAIL 

DEGREE  110 

1345 

203 

2592 

DVEC  3 

76 

1206 

143 

2308 

ADJCEN 

46 

881 

93 

1731 

LCA 

31 

648 

65 

1281 

BRANCH 

FACTOR  BY 

LEVEL 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RANDOM 

NOT  CALCULATED 

DEGREE  1 

1.15 

1.34 

1.57 

1.55 

1.79 

1  .93 

1.73 

1.59 

1.64 

1.71 

DVEC  3  1.07 

1  .26 

1.56 

1.60 

1.66 

1.70 

1.63 

1  .59 

1.93 

1.90 

ADJCEN  1.00 

1.00 

1.00 

1.14 

1.68 

2.00 

2.11 

2.37 

2.25 

2.52 

LCA  1.00 

1.00 

1.00 

1.14 

1.65 

1  .85 

1  .92 

2.15 

2.27 

2.32 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

RANDOM 

NOT  CALCULATED 

DEGREE  1 

1.89 

2.02 

2.01 

2.38 

2.41 

2.54 

2.66 

2.79 

3.09 

0.00 

DVEC  3  2.13 

2.26 

2.15 

2.27 

2.45 

2.45 

2.66 

2.77 

3.03 

0.00 

ADJCEN  i 

2.39 

2.33 

2.30 

2.20 

2.40 

2.55 

2.61 

2.86 

3.11 

0.00 

LCA  , 

2.39 

2.36 

2.44 

2.33 

2.42 

2.55 

2.87 

3.09 

3.57 

0.00 
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BRANCH  r*C7oR  BY  LEVEL 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RANDOM 

NOT 

CALCULATED 

DEGREE 

t  .14 

1.30 

1.47 

1.57 

1.72 

1.81 

1.85 

2.07 

2.30 

2.25 

DVEC  3 

1 .04 

1.17 

1.34 

1.59 

1.98 

1  ,9& 

2.04 

2.21 

2.22 

2.45 

ADJCEN 

1.00 

1.00 

1.04 

1.49 

2.02 

2.43 

2.39 

2.41 

2.25 

2.41 

LCA 

1.00 

1.00 

1.04 

1.47 

1.95 

2.21 

2.21 

2.29 

2.37 

2.50 

1  1 

12 

13 

14 

15 

RANDOM 

NOT 

CALCULATED 

DEGREE 

2.32 

2.57 

2.71 

3.14 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

OVEC  3 

2.49 

2.55 

2.70 

3.16 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

ADJCEN 

2.53 

2.63 

2.87 

3.22 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

LCA 

2.54 

2.70 

2.95 

3.38 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

N  =  15 

,  EDGES  = 

73(701) 

ORDER 

BRANCHES 

NODES 

avg 

max 

avg  max 

RANDOM 

NOT 

ENOUGH  SPACE  AVAIL 

DEGREE 

24 

75 

44  128 

DVEC  3 

21 

56 

38  96 

ADJCEN 

15 

32 

31  66 

LCA 

14 

28 

29  55 

BRANCH  FACTOR  BY 

LEVEL 

1 

2 

3  4 

5 

6 

7 

8 

9 

10 

RANDOM 

NOT 

CALCULATED 

DEGREE 

1.02 

1.07 

1.15  1.31 

1.42 

1.60 

1.71 

1.78 

1.98 

2.20 

DVEC  3 

1,02 

1.02 

1.06  1.24 

1.36 

1.60 

1.77 

2.01 

2.05 

2.33 

ADJCEN 

1.00 

1.00 

1.00  1.00 

1.08 

1.49 

2.03 

2.31 

2.40 

2.49 

LCA 

1.00 

1.00 

1.00  1.00 

1.08 

1.47 

1.93 

2.14 

2.14 

2.35 

1 1 

12 

13  14 

15 

RANDOM 

NOT 

CALCULATED 

DEGREE 

2.21 

2.52 

2.70  3.18 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

OVEC  3 

2.27 

2.46 

2.68  3.15 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

ADJCEN 

2.45 

2.54 

2.85  3.14 

0.00 

0.00 

o.oo 

0.00 

0.00 

0.00 

LCA 

2.49 

2.63 

2.91  3.44 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

N  s  30,  EDGES  =  57C30X) 

ORDER  BRANCHES  NODES 

avg  max  avg  max 


105 


N  s  20,  EDGES  =  1 33 ( 70X) 


ORDER  BRANCHES  NODES 
avg  max  avg  max 


RANDOM 

NOT 

ENOUGH  SPACE  AVAIL 

DEGREE 

52 

182  91 

329 

DVEC  3 

36 

119  67 

208 

ADJCEN 

21 

ai  a3 

76 

LC  A 

19 

35  ai 

66 

BRANCH 

FACTOR  BY 

LEVEL 

1 

2  3 

a 

5 

6 

7 

8 

9 

10 

RANDOM 

NOT 

CALCULATED 

DEGREE 

1.02 

1.03  1.07 

1.17 

l  .as 

l.as 

1  .67 

1.6a 

1.88 

1.87 

DVEC  3 

1.00 

1.01  1.05 

1.07 

1.25 

1.56 

1.57 

1.75 

1  .82 

1  .52 

ADJCEN  1 

.00  1 

.00  1.00  1 

.00  1 

.00  1 

1.01  1 

.13  1 

1.63  1 

1.97  2 

.30 

LC  A 

1.00 

1.00  1.00 

1.00 

1.00 

1.01 

1.12 

1.63 

1.88 

2.oa 

RANDOM 

NOT 

CALCULATED 

DEGREE 

1.80 

1.67 

2.17 

2.29 

2.52 

2.6a 

2.81 

3.10 

3.32 

0.00 

DVEC  3 

1.98 

2.11 

2.29 

2.a3 

2.58 

2.72 

2.8a 

3.07 

3.5a 

0.00 

ADJCEN 

2.a2 

2.52 

2.6a 

2.80 

2.65 

2.76 

2.99 

3.13 

3.a7 

0.00 

LC  A 

2.23 

2.ai 

2.53 

2.65 

2.77 

2.98 

3.19 

3.37 

3.91 

0.00 

N  =  25,  EDGES  =  90 ( 3QX ) 

ORDER  BRANCHES  NODES 

avg  max  avg  max 

RANDOM  NOT  ENOUGH  SPACE  AVAIL 
DEGREE  NOT  ENOUGH  SPACE  AVAIL 
DVEC  3  126  1158  251  2288 

ADJCEN  63  ail  132  819 

LC A  33  115  73  230 

BRANCH  FACTOR  BY  LEVEL 

1  2  3  a  5  6  7  8  910 

RANDOM  NOT  CALCULATED 

DEGREE  NOT  CALCULATED 

DVEC  3  1.20  1.35  1.58  1.85  1,80  1.93  1.98  2.05  1.79  1.8a 

ADJCEN  1.00  1.00  1.32  2.06  2.12  2.05  2.20  2.09  2.15  2.07 

LC A  1.00  1.00  1.32  1.97  1.91  2.08  1.98  2.03  2,08  2.08 

11  12  13  ia  15  1 b  17  18  19  20 

RANDOM  NOT  CALCULATED 
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DEGREE  NOT  CALCULATED 
DVEC  3  1.77  1.85  3.00  2.07 

ADJCEN  2.09  2.03  2.09  2.15 

LCA  2.02  2.03  2.09  2.22 


2.05  2.15  2.32  2.47  2.37  2.51 
2.26  2.16  2.32  2.3b  2.S2  2.53 
2.35  2.42  2.48  2.54  2.bl  2.75 


21  22  23  24  25 

RANDOM  NOT  CALCULATED 

DEGREE  NOT  CALCULATED 

DVEC  3  2.70  2.78  2.85  3.47  0.00 
ADJCEN  2.79  2.80  2.99  3.38  0.00 
LCA  2.89  3.08  3.37  3.83  0.00 


N  =  25,  EDGES  =  150(50X1 

ORDER  BRANCHES  NODES 
avg  max  avg  max 


RANDOM  NOT  ENOUGH  SPACE  AVAIL 
DEGREE  NOT  ENOUGH  SPACE  AVAIL 
DVEC  3  178  1328  340  2568 

ADJCEN  97  932  193  1829 

LCA  38  139  76  271 

BRANCH  FACTOR  BY  LEVEL 


12  3  4 

RANDOM  NOT  CALCULATED 

DEGREE  NOT  CALCULATED 

DVEC  3  1.09  1.24  1.29  1.U9 
ADJCEN  1.00  1.00  1.00  1.12 
LCA  1.00  1.00  1.00  1.12 


11  12  13  14 

RANDOM  NOT  CALCULATED 

DEGREE  NOT  CALCULATED 

DVEC  3  1.89  1.82  1.87  1.80 

AOJCEN  2.07  1.98  2.02  1.92 

LCA  2.08  2.01  1 .95  1.95 


5  6  7  8  9  10 


1.76  1.77  1.73  1.83  1.87  1.85 
1.40  1  .93  2.18  2.26  2.24  2. l  i 
1.34  1.86  1.91  2.08  2.01  1.99 


15  1 6  17  18  19  20 


1.86  1.92  1,94  2.03  2.08  2.16 
2.04  1.97  2.06  2.01  2.12  2.25 
1.92  2.04  2.04  2.20  2.26  2.37 


21  22  23  24  25 

RANDOM  NOT  CALCULATED 

DEGREE  NOT  CALCULATED 

DVEC  3  2.35  2.59  2.66  3.05  0.00 
ADJCEN  2.33  2.38  2.59  3.08  0.00 
LCA  2.59  2.65  2,94  3.29  0,00 


r 

- 1 

f  N  =  25, 

EDGES  =  210(70X5 

1  ORDER 

BRANCHES  NODES 

5 

avg  max  avg  max 

RANDOM 

NOT  ENOUGH  SPACE  AVAIL 

DEGREE 

NOT  ENOUGH  SPACE  AVAIL 

DVEC  3 

NOT  ENOUGH  SPACE  AVAIL 

.  ADJCEN 

87  723  169  1431 

LCA 

34  122  64  225 

BRANCH  FACTOR  BY 

LEVEL 

1  2  3  4  5 

6  7 

8 

9 

10 

RANDOM 

NOT  CALCULATED 

DEGREE 

NOT  CALCULATED 

DVEC  3 

NOT  CALCULATED 

ADJCEN 

1.00  1.00  1.00  1.00  1.00 

1.01  1.15 

1.48 

1.90 

2.06 

LCA 

1.00  1.00  1.00  1.00  1.00 

1.01  1.15 

1.43 

1.78 

1.88 

11  12  13  14  15 

16  17 

18 

19 

20 

| 

RANDOM 

NOT  CALCULATED 

DEGREE 

NOT  CALCULATED 

i 

DVEC  3 

NOT  CALCULATED 

r 

ADJCEN 

2.22  2.34  2.22  2.16  2.14 

2.13  2.14 

2.13 

2.05 

2.08 

t 

LCA 

1.94  2.05  2.12  1.95  2.05 

2.02  2.08 

2.04 

2.10 

2.27 

; 

21  22  23  24  25 

RANDOM 

NOT  CALCULATED 

DEGREE 

NOT  CALCULATED 

DVEC  3 

NOT  CALCULATED 

ADJCEN 

2.10  2.44  2.59  2.85  0.00 

LCA 

2.36  2.58  2.84  3.41  0.00 

1 

N  =  30, 

EDGES  s  130(30%) 

! 

ORDER 

BRANCHES  NODES 

avg  max  avg  max 

! 

RANDOM 

NOT  ENOUGH  SPACE  AVAIL 

OEGREE 

198  1504  389  2922 

, 

DVEC  3 

88  412  179  815 

ADJCEN 

63  691  134  1391 

* 

LCA 

31  57  71  120 

BRANCH  FACTOR  BY 

LEVEL 

1  2  3  4  5 

6  7 

8 

9 

10 

109 

kL__. 

- . 

RANDOM  NOT  CALCULATED 

DEGREE  1.36  1.73  1.97  1.82  1.89  1.97  1.92  1.86  1.9tt  1.93 

OVEC  3  1.21  1.48  1.64  1.74  1.94  1.91  1.92  1.95  1.92  2.03 

ADJCEN  1.00  1.00  1.30  1.86  2.16  2.21  2.26  2.22  2.32  2.25 

LCA  1.00  1.00  1.30  1.83  1.94  2.04  2.10  2.13  2.13  2.1o 


11  12  13  14  15  16  17  18  19  20 

RANDOM  NOT  CALCULATED 

DEGREE  2.09  2.18  2.08  1.97  2.05  2.04  2.12  2.15  2.37  2.32 

DVEC  3  1.95  2.05  2.16  2.33  2.25  2.22  2.24  2.36  2.30  2.34 

ADJCEN  2.28  2.21  2.25  2.20  2.31  2.32  2.21  2.16  2.31  2.39 

LCA  2.19  2.14  2.23  2.19  2.21  2.21  2.24  2.32  2.36  2.38 


21  22  23  24  25  26  27  28  29  30 

RANDOM  NOT  CALCULATED 

DEGREE  2.26  2.39  2.43  2.42  2.51  2.42  2.59  2.70  2.80  0.00 

DVEC  3  2.41  2.32  2.54  2.53  2.62  2.54  2.59  2.63  2.84  0.00 

ADJCEN  2.41  2.42  2.4 2  2.58  2.51  2.50  2.55  2.79  2.87  0.00 

LCA  2.46  2.51  2.52  2.60  2.72  2.81  3.03  3.12  3.41  0.00 


N  s  30,  EDGES  s  217(50%) 

ORDER  BRANCHES  NODES 

avg  max  avg  max 

RANDOM  NOT  ENOUGH  SPACE  AVAIL 
DEGREE  NOT  ENOUGH  SPACE  AVAIL 
DVEC  3  289  1825  555  3652 

ADJCEN  79  1995  163  3959 

LCA  32  76  70  148 

BRANCH  FACTOR  BY  LEVEL 

123456789  10 

RANDOM  NOT  CALCULATED 

DEGREE  NOT  CALCULATED 

DVEC  3  1.08  1.12  1.33  1.59  1.70  1.91  1.89  2.13  1.99  1.82 

ADJCEN  1.00  1.00  1.00  1.01  1.13  1.64  2.01  2.20  2.25  2.33 

LCA  1.00  1.00  1.00  1.01  1.12  1.63  1.88  2.00  2.01  2.08 


11  12  13  14  15  16  17  18  19  20 

RANDOM  NOT  CALCULATED 

DEGREE  NOT  CALCULATED 

DVEC  3  1.82  1.80  1.76  1.94  1.89  2.04  2.12  2.18  2.17  2.23 
AOJCEN  2.30  2.29  2.33  2.27  2.28  2.32  2.34  2.13  2.28  2.39 


HO 


LC  A 

2.18 

2. OR  2.13  2.20 

2.15 

2.18 

2.29 

2.24 

2.34 

2.40 

21 

22  23  24 

25 

26 

27 

28 

29 

30 

RANDOM 

NOT 

CALCULATED 

DEGREE 

NOT 

CALCULATED 

OVEC  5 

2.35 

2.43  2.44  2.51 

2.54 

2.62 

2.83 

2.94 

3.17 

0.00 

ADJCEN 

2.40 

2.37  2.43  2.45 

2.58 

2.62 

2.78 

2.96 

3.25 

0.00 

LCA 

2.50 

2.55  2.66  2.64 

2.80 

2.92 

3.06 

3.31 

3.75 

0.00 

N  =  50,  EDGES  a  304(70X) 


ORDER 

BRANCHES 

NODES 

avg 

max 

avg  max 

RANOOM 

NOT 

ENOUGH 

SPACE  AVAIL 

DEGREE 

NOT 

ENOUGH 

SPACE  AVAIL 

DVEC  3 

215 

1519 

394  2896 

ADJCEN 

46 

1003 

96  2003 

LCA 

30 

67 

65  132 

BRANCH  FACTOR  BY  LEVEL 


12  3  4 

5 

6 

7 

8 

9 

10 

RANOOM 

NOT  CALCULATED 

DEGREE 

NOT  CALCULATED 

DVEC  3 

1.00  1.00  1.04  1.05 

1.08 

1.22 

1.36 

1.34 

1.57 

1.65 

ADJCEN 

1.00  1.00  1.00  1.00 

1.00 

1.00 

1.00 

1.00 

1.05 

1.28 

LCA 

1.00  1.00  1.00  1.00 

1.00 

1.00 

1.00 

1.00 

1.05 

1.28 

11  12  13  14 

15 

16 

17 

18 

19 

20 

RANDOM 

NOT  CALCULATED 

DEGREE 

NOT  CALCULATED 

OVEC  3 

1.97  2.03  2.04  1.99 

1.83 

1.72 

1.68 

2.07 

2.16 

2.24 

ADJCEN 

1.67  1.97  2.19  2.35 

2.48 

2.52 

2.52 

2.45 

2.58 

2.55 

LCA 

1.66  1.88  1.96  2.09 

2.10 

2.26 

2.34 

2.34 

2.45 

2.45 

21  22  23  24 

25 

26 

27 

28 

29 

30 

RANDOM 

NOT  CALCULATED 

DEGREE 

NOT  CALCULATED 

OVEC  3 

2.33  2.49  2.64  2.70 

2.84 

2.98 

3.15 

3.34 

3.75 

0.00 

ADJCEN 

2.60  2.69  2.89  2.81 

2.93 

3.04 

3.18 

3.39 

3.83 

0.00 

LCA 

2.68  2.70  2.86  2.98 

3.13 

3.28 

3.48 

3.75 

4.23 

0.00 

OUTPUT  FROM  GRAPH  DIVIDING  ROUTINE 


1  1  1 


N 


EDGES  AVG  «  SUBS 


6 


8 


10 


12 


14 


AVG  SIZE  LARGEST  SUB 


n 

edq*s 

b 

3.4 

3.6 

3.6 

7 

2.5 

4.4 

5.2 

S 

1.7 

5.3 

7.1 

4 

1.3 

5.7 

8.7 

10 

1.1 

5.4 

4.4 

11 

1.0 

6.0 

11.0 

12 

1.0 

6.0 

12.0 

13 

1.0 

6.0 

13.0 

8 

5.4 

3.5 

3.5 

10 

2.6 

6.1 

7.6 

12 

1.4 

7.5 

11.3 

14 

1  .  1 

7.4 

13.4 

16 

1.0 

8.0 

15.4 

18 

1.0 

8.0 

18.0 

20 

1.0 

8.0 

20.0 

22 

1.0 

8.0 

22.0 

10 

7.0 

4.0 

4.0 

12 

4.2 

6.3 

7.7 

14 

2.7 

8.2 

12.0 

16 

1.4 

4.0 

14.8 

18 

1.4 

4.5 

17.3 

20 

1.2 

4.7 

14.5 

22 

1.1 

4.4 

21.8 

24 

1.0 

10.0 

24.0 

26 

1.0 

10.0 

26.0 

28 

1.0 

10.0 

28.0 

12 

8.8 

4.2 

4.2 

15 

4.3 

8.1 

10.4 

18 

2.5 

10.2 

15.5 

21 

1.7 

1  1.3 

20.0 

24 

1.4 

1 1 .6 

23.3 

27 

1.1 

11.4 

26.8 

30 

1.0 

12.0 

30.0 

33 

1.0 

12.0 

33.0 

36 

1.0 

12.0 

36.0 

34 

1.0 

12.0 

34.0 

14 

10.4 

4.6 

4.6 

17 

6.0 

8.1 

10.3 

20 

3.5 

lt.t 

16.4 

23 

2.4 

12.5 

21.1 

26 

1.8 

13. 1 

24.8 

113 


24 


26 


28 


27 

8.9 

12.9 

16.8 

32 

5.5 

16.5 

24.7 

37 

3.5 

18.9 

32.6 

42 

2.3 

20.5 

40.0 

47 

1.7 

21.1 

45.5 

52 

1.5 

21.4 

50.6 

57 

1.3 

21.7 

56.1 

62 

1.2 

21.8 

61.5 

67 

1.2 

21.8 

66.4 

72 

1.1 

21 .9 

71.8 

77 

1.0 

22.0 

77.0 

82 

1.0 

22.0 

82.0 

87 

1.0 

22.0 

87.0 

92 

1.0 

22.0 

92.0 

24 

19.8 

5.2 

5.2 

30 

8.4 

15.2 

19.7 

36 

4.5 

19.7 

30.1 

42 

2.6 

22.1 

39.3 

48 

1.7 

23.1 

46.4 

54 

1.3 

23.7 

53.4 

60 

1.2 

23.8 

59.8 

66 

1.1 

23.9 

05.8 

72 

1.0 

24.0 

72.0 

78 

1.0 

24.0 

78.0 

84 

1.0 

24.0 

84.0 

90 

1.0 

24.0 

90.0 

26 

21.8 

5.2 

5.2 

32 

10.5 

15.1 

19.8 

38 

6.1 

19.9 

30.0 

44 

3.9 

22.3 

38.5 

50 

2.7 

23.8 

46.5 

56 

2.1 

24.7 

53.6 

62 

1.6 

25.2 

60.6 

68 

1.3 

25.6 

67.1 

74 

1.2 

25.7 

73.0 

80 

1.1 

25.9 

79.7 

86 

1.0 

26.0 

85.9 

92 

1.0 

25.9 

91  .6 

98 

1.0 

26.0 

98.0 

104 

1.0 

26.0 

104.0 

110 

1.0 

26.0 

1 10. 0 

116 

1.0 

26.0 

1  16.0 

28 

23.6 

5.4 

5.4 

35 

10.5 

17.0 

22.4 

42 

5.5 

22.8 

35.5 

114 


? 

49 

3.4 

25.3 

'I 

45.5 

56 

2.2 

26.6 

54.1 

6  J 

1.6 

27.2 

ol  .3 

70 

1.3 

27.6 

69.1 

77 

1.2 

27.7 

76.3 

84 

1.1 

27.9 

83.8 

91 

1.1 

27.9 

90.9 

i 

98 

1.0 

28.0 

98.0  1  i 

105 

1.0 

28.0 

104.8 

112 

1.0 

28.0 

112.0 

119 

1.0 

28.0 

119.0  i  ■ 

JO  ] 

i 

30 

25.5 

5.6 

5.6  1 

j 

37 

11.9 

17.6 

22.7 

44 

6.7 

23.5 

36.0 

i 

51 

4.4 

25.7 

44.8  p 

58 

2.7 

28.1 

55.6  |! 

65 

2.0 

28.7 

63.0  fi 

72 

1.5 

29.3 

70.3 

79 

1.5 

29.5 

78.2  1 

86 

1.3 

29.7 

85.5  * 

93 

1.1 

29.9 

92.8  T 

100 

1.0 

29.9 

99.7 

107 

1.0 

30.0 

107.0  r 

114 

1.0 

30.0 

114.0  j 

121 

1.0 

30.0 

121.0  |; 

128 

1.0 

30.0 

128.0  ' 

135 

1.0 

30.0 

135.0 

40  il 

40 

34.8 

6.2 

6.2 

60 

7.9 

31.7 

49.0 

80 

2.7 

37.8 

75.9 

100 

1.5 

39.4 

98.6 

120 

1 . 1 

39.9 

119.8 

140 

1.0 

40.0 

140.0  1 

160 

1.0 

40.0 

160.0  ! 

180 

1.0 

40.0 

180.0  [ 

50  I 

50 

44.5 

6.5 

6.5  j 

75 

10.6 

38.7 

61.2 

100 

4.0 

46.7 

95.6 

125 

2.0 

48.5 

120.8  1 

150 

1.4 

49.5 

148.7 

175 

1.1 

49.9 

174.9 

200 

1.0 

50.0 

200.0 

225 

1.0 

50.0 

225.0 

250 

1.0 

50.0 

250.0 

60  1 

115 

400 

1.1 

99.9 

399.2 

450 

1.0 

100.0 

450.0 

500 

1.0 

100.0 

500.0 

550 

1.0. 

100.0 

550.0 

*  *  * 

*  * 

*  *  * 

****** 

OUTPUT  FROM 

************* 

DIVIDE  AND  COLOR  PROGRAM 

*  *  * 

*  ★ 

*  *  * 

****** 

************* 

N  s  25, 

EDGES  = 

6(  2%) 

ORDER 

BRANCHES 
avg  max 

NOOES 
avg  max 

RANDOM 

6 

6 

12  12 

DEGREE 

6 

6 

12  12 

DEG  3 

6 

6 

12  12 

AOJCEN 

6 

6 

12  12 

LCA 

6 

6 

12  12 

N  =  25, 

EDGES  = 

12C 

4%  ) 

ORDER 

BRANCHES 

NODES 

avg 

max 

avg 

max 

RANDOM 

12 

12 

23 

24 

DEGREE 

12 

12 

23 

24 

DEG  3 

12 

12 

23 

24 

ADJCEN 

12 

12 

23 

24 

LCA 

12 

12 

23 

24 

N  s  25, 

EDGES  = 

1 0  C  6X) 

ORDER 

BRANCHES 

NODES 

avg 

max 

avg  max 

RANDOM 

17 

IS 

34  36 

DEGREE 

17 

19 

34  36 

DEG  3 

17 

18 

34  36 

ADJCEN 

17 

18 

34  36 

LCA 

17 

18 

34  36 

t 


I 


N  =  25,  EDGES  =  24(  8X) 


ORDER  BRANCHES  NODES 


avg  nu 
RANDOM  21  25 


avg  man 
a5  55 
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DEGREE 

22 

46 

46 

90 

DEG  3 

21 

27 

44 

55 

ADJCEN 

21 

23 

43 

54 

LCA 

21 

23 

44 

54 

N  s  25, 

EDGES  = 

30(10%) 

ORDER 

BRANCHES 

NODES 

avg 

max 

avg 

max 

RANDOM 

25 

45 

53 

91 

DEGREE 

27 

52 

57 

107 

DEG  3 

2a 

31 

52 

68 

ADJCEN 

23 

30 

52 

67 

LCA 

23 

30 

52 

68 

N  =  25,  EDGES  =  36(12*) 


ORDER 

BRANCHES 

NODES 

avg 

max 

avg 

max 

RANDOM 

51 

945 

108 

1893 

DEGREE 

34 

191 

72 

380 

DEG  3 

26 

94 

59 

189 

ADJCEN 

24 

35 

57 

82 

LCA 

24 

35 

57 

82 

N  =  25, 

EDGES  = 

42(14%) 

ORDER 

BRANCHES 

NODES 

avg 

max 

avg 

max 

RANDOM 

84 

1012 

172 

2021 

DEGREE 

43 

208 

91 

425 

DEG  3 

29 

86 

64 

177 

ADJCEN 

29 

170 

65 

349 

LCA 

25 

54 

58 
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N  s  20,  EDGES  =  38(20%) 


order  branch®* 
avg  max 


nodes 
avg  max 


Logarithm  of  max 
branches  nodes 
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DEGREE 

46 

363 

96 

731 

5.89 

6.59 

DVEC  3 

27 

137 

60 

278 

4.92 

5.63 

ADJCEN 

23 

82 

52 

166 

4.41 

5.11 

LCA 

20 

33 

46 

70 

3.50 

4.25 

N  =  20, 

EDGES  = 

47(25X3 

order 

branches 

nodes 

Logar i t hm  of  max 

avg 

max 

avg 

max 

branches  nodes 

DEGREE 

45 

127 

94 

248 

4.84  5.51 

DVEC  3 

28 

93 

62 

182 

4.53  5.20 

ADJCEN 

21 

43 

52 

87 

3.76  4.47 

LCA 

20 

27 

48 

63 

3.30  4.14 

N  =  20, 

EDGES  s 

57  ( 30% ) 

order 

branches 

nodes 

Logar i t  hm 

of  max 

avg 

max 

avg 

max 

branches 

nodes 

DEGREE 

49 

187 

99 

380 

5.23 

5.94 

DVEC  3 

32 

164 

67 

320 

5.10 

5.77 

ADJCEN 

25 

109 

56 

223 

4.69 

5.41 

LCA 

21 

55 

48 

117 

4.01 

4.76 

N  s  20, 

EDGES  s 

66(35%) 

order 

branches 

nodes 

Logar! thm  of  max 

avg 

max 

avg  max 

branches  nodes 

DEGREE 

112 

1119 

219  2198 

7.02  7.70 

DVEC  3 

52 

308 

104  608 

5.73  6.41 

ADJCEN 

38 

513 

80  1015 

6.24  6.92 

LCA 

23 

62 

51  123 

4,13  4.81 

N  s  20, 

EDGES  = 

76(40%) 

order 

branches 

nodes 

Logari thm  of  max 

avg 

max 

avg 

max 

branches  nodes 

DEGREE 

NOT 

ENOUGH  SPACE  AVAIL 

DVEC  3 

59 

267 

116 

524 

5.59  6.26 

ADJCEN 

30 

292 

65 

593 

5.68  6.39 

LCA 

23 

285 

52 

568 

5.65  6.34 

N  a  20, 

EDGES  = 

85(45%) 
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order 

branches 

nodes 

Logar i  t  hm  of  max 

a  vg 

max 

avg 

max 

branches 

nodes 

DEGREE 

92 

635 

173 

1206 

6.45 

7.10 

DVEC  3 

5o 

425 

109 

818 

6.05 

6.71 

ADJCEN 

37 

379 

77 

749 

5.94 

6.62 

LCA 

23 

64 

49 

128 

4.16 

4.85 

N  =  20,  EDGES  =  95(50%) 


order 

branches 

nodes 

Logar ( t  hm  of  max 

avg 

max 

avg 

max 

branches 

nodes 

DEGREE 

101 

811 

188 

1608 

6.70 

7.38 

DVEC  3 

57 

367 

107 

676 

5.91 

6.52 

ADJCEN 

34 

501 

71 

980 

6.22 

6.89 

LCA 

28 

584 

60 

1150 

6.37 

7.05 

N  =  20,  EDGES  =  104(55%) 


order 

branches 

nodes 

Logari thm  of  max 

avg 

max 

avg 

max 

branches 

nodes 

DEGREE 

106 

431 

194 

777 

6.07 

6.66 

DVEC  3 

69 

440 

131 

844 

6.09 

6.74 

ADJCEN 

34 

258 

70 

497 

5.55 

6.21 

LCA 

23 

79 

48 

156 

4.37 

5.05 

N  =  20,  EDGES  s  114(60%) 


order 

branches 

nodes 

Loga r i t  hm  of  max 

avg 

max 

avg 

max 

branches  nodes 

DEGREE 

99 

9b9 

iao 

1696 

6.88  7.44 

DVEC  3 

62 

576 

114 

1029 

6.36  6.94 

ADJCEN 

37 

561 

75 

1088 

6.33  6.99 

LCA 

27 

438 

56 

855 

6.08  6.75 

N  s  20, 

EDGES  s 

123(65%) 

order 

branches 

nodes 

Loga rithm  of  max 

avg 

max 

avg 

ma  x 

branches  nodes 

DEGREE 

86 

765 

155 

1351 

6.64  7.21 

DVEC  3 

55 

466 

101 

832 

6.15  6.72 

ADJCEN 

28 

188 

57 

365 

5.24  5.90 

LCA 

21 

53 

43 

103 

3.97  4.63 

N  s  20,  EDGES  s  133(70%) 
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order 

branches 

nodes 

Logar i t  hm  of  ma 

avg 

max 

avg 

max 

branches  nodes 

DEGREE 

74 

419 

130 

748 

6.04  6.62 

DVEC  5 

47 

282 

8b 

498 

5.64  6.21 

ADJCEN 

21 

57 

44 

114 

4.04  4.74 

LCA 

20 

33 

41 

62 

3.50  4.13 

N  *  20, 

EDGES  = 

142(75*) 

order 

branches 

nodes 

Logar i thm  of  max 

avg 

max 

avg 

max 

Drenches 

nodes 

DEGREE 

60 

295 

103 

499 

5.69 

6.21 

DVEC  3 

38 

226 

69 

403 

5.42 

6.00 

ADJCEN 

22 

66 

43 

118 

4.19 

4.77 

LCA 

19 

37 

39 

71 

3.61 

4.26 

N  s  20 

,  EDGES  = 

152(80*) 

order 

branches 

nodes 

Logarfthm  of  max 

avg 

max 

avg 

max 

branches  nodes 

DEGREE 

45 

143 

76 

237 

a. 96  5.47 

DVEC  3 

30 

71 

54 

118 

4.26  4.77 

ADJCEN 

20 

46 

39 

81 

3.83  4.39 

LCA 

19 

27 

37 

49 

3.30  3.89 
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